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ABSTRACT 


A tedmlqiw Is dsvsiUvad for the estiaation of total daily insolation on the 
tsaais of data derlvabls fh>a qperatioaaIpclai>ortdtingaBtfll]ita& Althou^ aurfhca 
insolation and aataorologLcal observations are used in the devali^Baiit, the 
algoritiiB is oonatrained in ig^oatlon the inft'aquent daytiae polai>«rtitflr 
coverage. 

A lineer regraerion ^proadi is ueed to relate o b ser ve d Insolation to the 
cgjtloal tJTloknafH, e x pr aese d water ve^ior slant patti and effective cloud 

trensBittanoe (p^qoMog cloud aaouint and type). Reip'eaBlon analyses relating 
these peraBetere to the total daily in'cIa'.«Qc result in daily estlBates for each 
hour. EstlBatas o oneap ond ing to the u.Bes of satellite psaaage are weighted, 
according to their proxlaity to looeiL solar noon, and avera^ to provide a aii^ 
daOy iiBolatlexi estlaate. 

Prior to final tuning of the ooeffkdiKits, an interactive test of the proeedire 
was conducted on inc^p^enb data, with GOES satellite iBagery as the data 
soiace for polarMXtiter tiBes. Results were suooesBf\iL in ter bs of the ratio of the 
standard deviation of the redduala to the ( S we rv ed Bean Class than 0.2 ttreehold) 
and conation coefndaits ncceeding 0.80. After tuning the doud 
paraBeterizatlons and develi^dng a final set of regression coefficients, checks on 
in<tep«ida3t surface data revealed ttiat auidltlMial iaproveB«its can be Bade by 
introducing a Uaa reduction tedinlgue. 

An autoBated procedire for doixl type daastficat'iMi baaed on the 
Bultiapectral infrared sounder data was stipulated fTom oim data swath and tested 
oa another. This type of dcul ciajwtftoflWon ftom operatlcxial data shows 
considerable proBise for fixtise ^plication. 

This document has been prepared as a part of the FY81 Yield Model 
Development Project. 
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1*0 niTROPocTioii 

1.1 

Vm ApteuStiM nd R«wura« Ita fio ba t y Sirv^ ttirougb kmrotptm 
Rtaote Stariag (AgRTiTARS) pf^ogna iflan ctSkt tar •BMantwi of oeaar mdlift&oii 

tb§ Bxtkoa bMtd on oporatloniL wlMilltn data. SoadatlQn 

and otbir quMttttw ara naadad botti Par aarfy aaniing of dbui§m ^tadtiog 
roaouroa quautg Mid Par oonnodttgr produotlon Pa ra oa at a . To naat raqift<aBaDt% 
tadmiquaa Par aafcLnattan of Inaoiatton nuat ba datal^Mal Md than put Into 
oparatlonal ig)plllaatton Par uaar taaPng ovar alpiiflottit avdotdturaL im 
wortd. 

Pwnlatlon ia an AgRSTARS qunttty nora oiloaaljr aanodatad with Input 
to orop yie3d nodeila CPar oon modify production P ara oaat a ) Vtma tar othar 
ai^fllMtlona, but tha Inaolation alao ia a oonfaibutor to the aatLaatlon of othar 
quantitiao ( 04 ., aoQ. aoLatira) that raiata to aartr aanihig of otangaa of 
tdgnifloanoa to agrloidtura. (Uobal oovaraga raquiraattita can Par at laaat one 
tadnfyia that will provida looal aatinataa of <hdfy total inacaatlon Par any r^lon 
of thegloba. 

Moat of the r aoan t efPoit to extract inPomation on inacdation Pion 
aataUita data has o on o a idra ta d on the mai^ of radianoa aiaaairananbB of 
raflaotad acilar radiation obtained Proa gaoaynohronoua aata3Utaa. TlMaa pHatforaa 
ofPar P^acpiant (half-hoir) faigh-raaoSution aaaauraaanta of aoSar radiation Pafling 
to raadi tha aiaPaoa. Althafyh it ia poartMa to oopa with the data handling 
profaiLeaa Par routine prooaarii^ of high-renoaution digital data, not aU gaogmphloal 
araaa of intaraat ar% praaentfy oovarad. Data Pi^a tha ^aratioul poiar-orbitiiv 
aataUitaa provida an altamatlva with regular global oovaraga. Tha kay quaablon ia 
whether or not tha inPBquant obaervatiooa ftoB tha polaiMxtitar idUauPdoa Par a 
apaodfioation of the daily t^al insolation, whidi di^Mndi oritioally on tha 
olixidinaaB throu^ 1^ oantral auiilit hoixa. Utilization of tha aUltiapaotzal 
oa p a b Di t iee of tha polar cctitare night ooapanaata Par the inPaqunoy of 
obamratioaB at a glvm looatioo. 

It h oontaaplatad that tha Inaolation algoritba wflL initially ba 
iapSaaentad intaraotlvaly, in oonjonotion with praolpitation aatiaation tachnkpiaa. 
OptiauB <tata dMmotariatiQa Caxh u high readution, or noontiaa data) aay not 
ba availaUa. TharaPora, tha aatiaation tadmbpia nuat ba flaxihla «»ugb to aal» 
tha beat use of whMaver inPoraatlon ia availabla oparattonally. 
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Ottiwtivis 

Th« ai^ of thia atudly haa baan ttm tevaaopBant of a 

tadttiiqua fbr aafl aat to n of tha <tid]y tottL taaoaatton uilng date aoqidrad fhoa 
op«mtioQal polar-oittting a^anitaa. Aa a fSrat {iaa^ aa alforttha ftr tba 
ta^siLi^ b daatgnad on tea baoii of < ^ amr a tioni at tea aurteoa. Atepted 
paraaateilKatioaa ara to ba aadaftt>B ^cmtlotMa. aatamta tete - - poatfdy froa 
atitttopaotrAl ataoaplMrio aoundar date. IJbmmietnBy, tea paraaatera of tea 
algoritea ars aatlaatad by intaraotiva atelyaia of routina aataHtba laagwy of 
Bodarata naoaution, auppiaBantad tqr oo na a i n n al axfboa Utfcraatlon or prodiMte 
from tt» opaodtonal aataHlta aaaadar, HialntmdadteatteadaivadalcortteB ba 
aada avaflabLa for roirtiM taattng. 

1*3 Soopa 

Tfala r^wrt Innlwlaa a vary Uaitad dtootaalnn of badqpraund, tete 
aouroaa» and iraUainary anaiyaaa in tea aatabUahaant of tea ^ni^oadu A 
diacuBBi^Mi of tea ai^roaoh ravada llBitetionB teat hava a having on tea raauilte 
obteinad and on tea p roa p a ^ for f\itura iBprovaaante. 

RaeKdte teftar aodlflMtiona) of tea r^aaeion ana^aaa ara aua aartzad 
for aadi a a aa o n fh>a 2} ytmn of date ovar tea Unttad ^ataa. TabuOatlons of tea 
lataat waigiiUng ootfOoiante anaUa i^i3ioation of tea aateod* A briafaooountof 
indapandant date taata ia glvan, aa wan aa a dlaraaadm of tea faawlhnity abudy of 
tea diract aetiaation of ctloud pvanatara fix>a aultiapaotrai infterad aoundar 
data. 
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24) BACKQROOMD 

H«ebodi ttet tatvt b&m te ¥C k « d to provite iitlaatio of tht tdtA 
on tbo bMii of attiUta^tarlirif tnftraaittnQ htf« oapiiMlMd «ith«* w tapMool 
■Wwooh or > Piaifllflad phaatori. aodoL Rijorouo mdiaitiy traniftr ooaputatton% 
ovir tte oo ttr t solar apsotrua sad indudias all ord«a of soattartr^ am not 
praottosMs far routtos ^^piioat^ Pia1;b«*aors» adlpartin«itoptioad.fafarastion 
far ttm udattiig staosphart aoA odouda at u t teaiy tiaas Md looitiona art not 
avafla^ la any oase. botti aapMoal and ptayafoal aoMs ttat ara a^iiad 
tfSI aaka uaa of poraaatari ttMt ara lUtaly to ba baaad on QHaatoloBtoil 
Infaraation* Whaoavar poadhla, nowavw*, tha aigntfloant vartationo of olouds and 
watar vapor win ba daptotad on tba baafa of oponatlontl Infaraation. Cloud 
pr^artlao and total watar vapor, objaotlvaiy daaort.Dad, baooaa tha laportent 
ooaponants of any taotniqua. Raaota aanefag data aicht ba uaad to infar audi 
proparOas of Intaraat indiraotly, or a4hk ba n aa d ad to obtain aoaa rdatad 
Infaraation diraotly (a.g., doud raflaotanoa). 

E«4y affarts to aaka uaa of aataTMta data In aatiaation of ataifaoa Insolation 
(of. Fritz, at d.^) daaonatratad ralationahlp batwaan aataillta aaasuraaanta 
and ataaapborlo attanuation* TbadipdfloanMof doudaaounthasbaanavidantin 
an affarts (of. Hraaon at aL^, Qvdnn3, SQlto uid Vondar Haar^} to aatahllah 
pMvtgatartzad a xp raaaiooo. Anrasiiltaarodapandantonthaapacja-tlBaaoalao.and 
solar angina oonaldarad and, on tba variable atnoaphatio attanuation} thay an 
nnkad IndiraoUy to aurfaoa raflaotanoa. Not an raoant aodals uaa aataillta data 
but d^and rafaar on oonvandonal aurfaoa obaarvations and pmnatarizations (of. 
Atwatar and Brown^, Atwatar and Ban^). Thasa atudliis ravaSl tha laportanoa of 
doud typa Infaraation. Vlfa adi^itationa, aataillta data oould ba uaad far ir^ut 
Infaraation If doud typas oodd ba aatlafaotoiily ohnaoterizad. Another 
nuaadoal aodal that has baan ^pllad «ctaneivaly la that of Laoia and Hanaan^. 
Again, far tha diraot adaptation of aataTHta data, trandatioca to optioal 
poraaatarB ara naoaaaaiy. 

Currant i^proadias noat daamotartstic of tha uaa of sdar raflaotanoa 
BaaauraaantB faoa gaoarndinoooua aataOltaa a.a tha aapirloal rairvwfnn rodal 
(Tarplay^) and tha Maplifiad phyrioal aodatL (Gautiar^). Both adhods raqdra 
tsightnaas asMiaaaants far bofa Oloudf and oloud-faea badcfi'oundo. The Gautiar 
nodal raquiras apaoIfUation of ataoapbario attanuation ooefOoiaiito and otbud 
abfMTption ooaffloianta, whanaaa tha TanAay nodal ooefOdants «a baaad on laoat- 



The basic pu^xsae of sonitGring dcud rdlactance is to determirw the 
reduction in inaoHatlcxi. The more fhequenOy ttiis is done, the better wDl be the 
estimate of the daily insolation eeqpecially Mhen the observationB are made dose to 
the central part of ^ day. 

An dtematlve to meaeiring reOeotance is to infer (firom measurements) the 
transmittance of available radiation to the surfhce. This ai^roach is necessary 
whm insufficient daytime reflectance measurements are available, as is the case 
forfliisatudy. Ultlmateily, it is neceoary to account fbr transmittance in any case 
wlMn working with ane-edded satellite measiremaits of upweQling radiance. 
(AnoflMT dtematlve is to base the estimate of insolation on a m^mured response 
of the surface, through its tempmatire, but such a pro^mm of remote 
measurement requires information on cloudiness and surface diaracterisUcs as 
welL Furthermore, it is preferable instead to estimate the impact of insdatieai on 
the thermal respo ns e). 

Previous teduxiques have revealed the difficulties both of the statistical 
approadi to the data and of the digital handling of the satdllbe data. Some bias in 
results is anticipated for estimates of octremes. As for the development of tlw 
application of satellite data, it win require good quality ground truth data of the 
insdaticxi at locdities that are representative of areas covered by remote sensing 
measurements. 
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3.0 DATA 

3.1 SOLMET Data 


SOLMET. houdyScJar Radiation »SurfiMg Mete oroilDidcalObBervatiMM. 
is an aot±ived dataaet, on ai^notie tape p rod u ced bf the Environfflental 

Data and Difb»*mation Service of ttw Naticxal Ciimate C«itei^ SOLMET coiBtdnee 
houriy incoming solar radiation data with the corresponding meteorological 
c^jeervation over the 3d^tion NOAA Solar Radiatim Netwocic into one 
compr^ienaive dataset. Two and one-hailf yeen at data ftvm twenty of the 
SOLMET statiixs were used in tiito analyi^ Fl^pre 1 IDuatrates ttiese locations 
geog^I^caTly, whOa Table 1 provides pertinent station information. 

For the insolation, use was made of the meesired and edited gik:^ 
(total direct and diffuse) solar radiation on a horizontal surface. For tiie 
extraterretizial radiatiixi, listiz^ of the computed available solar &aergy received 
by a hori^xital surface at the t<^ of the atmo^^ieie (based on a solar constant of 
1377 J/m^ * s), were used. The data are Gagged to indicate source, if other than 
obsoTved irradiance, or to indicate corrections. The meteorological data used 
c(»)tisted of scrface preesire and dew point as wen as detailed cloud layer type and 
V amount informatiicxu (Visibility <^servatiot)s were not used in this study). Table 2 

sum marizes most of the SOLMET data actually incorporated into the analyses. 

3.2 Satellite Data 

The SOLMET data with the si;q>potting houiiy (and 3-4io(jriy) surfkce 
meteorological obeervations cover only very small portions of the area under 
routine satellJte siaveillance. Each polaiMrhiter can view almoet any ar^ twice a 
day (four times for tbs system pair of operatitxal polar orbiting sateOttes). Of 
particular interest are the TIROS-N, NOAA series of satellites'* ‘I. The operational 
satellites c<»itain two scanning systems with visible and infrared radiometric 
measuremoits: (1) the Advanced Very Hl^ Resoluticm Radiometer (A VHRR), and 
(2) the Hi^Resolution Diftared Radiometer Sounder (HIRS). The latter scanner is 
part of the TIROS Operatinxial Vertical Sounder (TOVS) system aiKl plays a major 
rde in the gaieration of TOVS retrieved products. 

The A VH R R data are available in Gve differmit spectral ranges, two of 
which are in tte short wavelength regk»i (0.58-0.68 urn and 0.725-1.10 urn) and 
three in the infrared region G.7 urn, 11 urn, 12 um). Spatial resolution near nadir is 
about 1.1 km. It was inteided that daytime short waveLmigth AVHRR data 


5 



FIGURE 1 NOAA SOLAR RADIATION NETWORK 
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TABLE 1 

SOLMET STAnONS USED IN STUD! 


STAHON # 

NAME 

LAT 

LON 

03937 

Lake Charies, LA 

30.07N 

93.13W 


ColuBhia, MO 

38.49N 

92.13W 

12839 

Miami, FLA 

25.49N 

80.17W 

12919 

^traavina, TX 

25.54N 

97^W 

13722 

Rakd^ NC 


78.47W 

13895 

Montgomery, ALA 


86^4W 

13897 

NaahvULe, TENN 


86.41 W 

13985 

Dodge City, KAN 


99.58W 


Caribou, ME 


68.01V 

14742 

Biriington, VT 

44^N 

73.09W 

14837 

Madiaon, WIS 

43.08N 

89.20H 

23023 



102.11W 


Bismarck, ND 

46.46N 

1 00.46 W 

24143 

Great Fa^ MON 

47^N 

111.22V 


Fairbanks, Alaska 

64.49N 

147.52V 

93734 

Steriing, VA 

38.59N 

77.28V 

93805 

Tallahassee, FL 

30^3N 

84.22V 

93819 

lhdianapolia,IND 

39.44N 

86.16V 

94823 

Pittidjurg, PA 

40.30N 

80.13V 

94918 

Omaha, NEB 

41.22N 

96.01V 


ELEV 

19 

277 

8 

12 

137 

68 

186 

795 

195 

112 

271 

872 

511 

1118 

143 

87 

18 

244 

371 

404 
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TABLE 2 

TYPES OF SOLMET DATA USED IN STUDY 

WBANSTAHON NUMBER 

YEAR 

MONTH 

DAY 

SOLAR HOUR 
SOLAR MINUTE 

ESTRATERRESHAL RADIAHON 
DIRECT RADIAnON A FLAG 
INDIRECT RADUTION & FLAG 
OBSERVED TOTAL INSOLAHON AND FLAG 
STAHON PRESSURE 
STATION DEW POINT 

SKY CONDITION (cLrar, acattered, broken, overoaat, or obacired) 
TOTAL SKY COVER 

LOWEST CLOUD LAYER AMOUNT AND TYPE* 

SECOND CLOUD LAYER AMOUNT AND TYPE* 

THIRD CLOUD LAYER AMOUNT AND TYPE* 

FOURTH CLOUD LAYER AMOUNT AND TYPE* 

TOTAL OPAQUE SKY COVER 


*Glven every three hours, baaed oa 16 genetic clixid types or obscuring phenomena. 
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at \3m tdg^ reaoSutioa would ba uaad only ae a maaaa of daaivdng and ^Motdng 
doiH apadCloatioiiB, and may ba avaSabla only onoa par day. For r^ular analysis 
it is mora likily that AVHRR data that ara storad at radiKsad resolution (Globdl 
Area Covarega), with a subpoint rasolution of 4 km for both the inftsrad and 
availsbls visible data, would ba used. The raduoad resolution, with only a third of 
the original aoanlines used, oonalats of 409 spots in met soahlina. After 
IdantifLoatim of AVHRR pixels within soma qpaolfled grid box, ttireeholding and a 
histogram analysis could be patfbrmad for cloud amount, bri^itneaB, and 
tamperatire. The data may be analyzed interactively fl?om multichannetL Imagery. 

HIRS data, with an optimum spatial resolution of about 17.5 km, 
re p rese nt a potential aouroe of data for routine analysis of cloudiness within the 
general areas of interest. Eadi aoan Una contains 56 scanspots of data for ea^ 
channeL Every 40th acahllne there usually is a 3^line data gap for calibratim. The 
HIRS instrument includes 19 Inftsred radiometric channels qpaced fh>m the near- 
infbarad window regim around 3*7 urn through the central portion of the 15 urn 
band of carbon dloxlda. A sin^ short wavelength channel Is omtend around 0.7 
micrometers. This channel is avaOsble, under suftlclmt suhli^it, as a general 
souroe of brightness data to supplement daytime inferences made fh>m the inftsred 
data. However, it is likely that AVHRR data would stilL be required for the initial 
development at doud amount ^^edflcati^is. If HIRS and AVHRR data are applied, 
say, to a 2.5 X 2.5 degree grid box with a subgrid of 0.5 degrees (25 subdividmB), 
then several HIRS datapoints can be averaged for each subdivision. Within each 
subdivision, it would be possible to obtain statistics of AVHRR data (over 150 
points archived at 4-km resolution}. 

AVHRR and HIRS data available an digital tape include lcx»tiixi and 
caUbratiEm informatiocu However, limb correction algcrlthms must be ai^plied for 
standardizaticsi prior to analysis (spectral varlatixxiB in infimred Slant path 
attenuation with zenith an^ are thereby accounted fcr). HIRS data swaths ftom 
successive orbits begin to overlap poleward of about 35 depees latitude. 
Equatorward of that latitude it is possible to have data gaps to tlM extent that a 
particular spot may not be oovered by the orbital swaths for that mode (ascending 
or descoidlng) on that day. Also, on successive days, the satellite subtrack shifts 
about 8 decrees of Ixxigitude ralatlve to a fixed equator croaslng. 
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With an additional reliction in apatial oovaraga, standard TOVS 
products are routixMily available far qpMifled retrieval boxes. For eadi box these 
products incQu(te dear^iObluon HIBS radiatioes, retrieved temperatures and 
retrieved mixing ratios appropriate to q>eci£led pressure levels. The total 
predpltable water may be Inferred i>iom the mixing ratios. 
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4.0 APPROACH 

4.1 Planned Analvaae and T<ats 

4.1.1 General CaMbrainta 

The apiroaoh axk^ited for ttie uae of data obtalnahle oparatLonally 
ttom the poQap^rbttiz^ NOAA aataUites has evolved on the preDlae that ttw 
preferred r^Mfettlve hi^x^aaolutlm vIMhLe rafleotanoe data wlIL not be avafiabila at 
^ptinuB tlaee the day. TUa lialtatlon of the avaUabile po]ar.crbiItar data 
doaa not altar tfaa major (^^Botlve of any aucoeeafUI ^proadi, naaaily, to aooount 
for the iapaot of dloud attenuation on the daily inaolation. Toward that end it ia 
BOflt d e si r a b l e to deal with a physical modal that aooounta for the radiative 
prooeeses. A rigorous ladiatlve transfer treatment would not be acceptable on 
operational gi'oundey both fi?om the point of view of available infcrmatifxi and 
compu t a t io n al demaikb. Greatly simplified physioal models can be introduced ae 
long as adequate perameteKizations are incorporated. Componenta of the model 
dK3uld be "obeervabLe" thorn the satellite platform. However, when reflectenoee 
themselves are not incorporated in the mo<^ then the link b^ween the remote- 
sensing data and the total tranamittance of solar nuiiation to tte lower boundary 
autfhoe must be achieved indirectly ttirough more parameterization. 

4.1.2 Regreaeion Analyala. Depmdent Data 

The di fflo Gty of aatiafactocily accounting for atmoeplMric 
transmittance thorn (xie^ided measurements dictated the ^rrfWaHnn of a re{p:^eaBion 
aptproach, using actual <^)eervations near the sirfaoe. Ihitially, the ind^p«ident 
variables in the i^egreesion approada are selected thorn surface data, but are rixseen 
to r^resant informatian that can be infen^ed thorn polai>orbiting data. 

The i^egreasion expreasLon, the other hand, is selected so aa to re p r e a en t (in a 
very gmeralized way) the {lyrical proceaaes. Success or failure at the regr^eesLon 
^proaoh depmds oi the suitability of the avaDable data and on tlM writahruty of 
the {diysicalrepresentativenesB of the regreasicwi expreasLcm. 

During the design stage at the regreaeion analysis (with depmdent 
data) observaticxB of total daily inaolatitx), along with oompubed insclatiixi inddoat 
on a iMrizontal surface at the tc^ of the atmo^iiere, are used to desribe the 
dependttit vatiablB. doiri obeervatdxxB (amount and type) and aurfaoe dew point 
temperatures an uaed to describe indeperxient ^^riabLae indicative of 


11 



trtmmtttanott ttrough doudy and c3oud-A:M atmoephariB (viVti wiabaa watar 
vi 4 )or). Exparianoa with thaaa data typaa alao provldaa ttia infioraation naadad for 
laproving tha modal or for altaring tha dataaat whan poaalMa. Tha propaa atarta 
wiO) axaolnation of hourly data* ainoa both tha d^andant and indapmdmt data 
aouroaa are avaflahla at that lh«;uanoy. Analysia of houdy data alao makaa It 
aosiartoaaaeaatteraiLaof diangaalnthaaoklarzantthan^. N«ct* tha total ddly 
Inaolatlon la axaalnad In ralation to apaoUlo houriy aouroaa of Ir^lapmdent data* 
since tha data win be acquired ody at aeveral qpadflo hours (Arom poiLBOMxtlting 
aataHitaa). 

4.1.3 Jhaaowadant Data Teats 

Onoa tha regr eaa l on analyaea with the dependant data are 
astablMied* axpreasiona must be tasted. One type of test la to make use of tha 
SOLMET data themselves for periods not previously Inducted in the design andyaia 
with d^Mtidant data. For that program, two monthe of data during the 
agriculturally signifloant spring period were reservisd for testing. 

After pro 0 ?eaBlng to the stage of inootroratlon of operational 
satellite data aa the aouroe of indepeiwtent data* the first testa are in terma of an 
interactive aneilyais of simple aatdllte Imagery, primarily fhom ain^ channd. 
infimred radiometalc meaaurem«its* as auppLemrated during the day by vistble 
refiectance meaairements. These maaau^ments are typical of operational 
soarmece aboard Mm polarvcrfaltlng aatSQites, althoudi the high resolution scanners 
<xi the most recent pdaivxtitera are not limited to a ainds visible and a aln^ 
infbared channeL Neveiihelaas* to simulate palaiMxtlter coverage* satellite 
Images were selected thorn fdlrdtec InftBrsd GOES data* almg with occasional 
visible Imagery, to rep r ea en t imagery that mi^it be available to the operational 
analyst for interactive scrutiny. Generally, only three images were selected to 
repi^asent the eariy morning-early evening paseaa and Mm mld-attemoon paaa of the 
NOAA aeries polar orbiterB. For this test, doud perameters had to be estimated 
for local areas immediatdy aroutKl specific pyranometer statiniB* whUs the 
molstixv psTameter was estimated fhom NMC hemiqherlc ( 6 -hour) reports on 
surface diarts. 

Ih addition to the interactive test of the me^Kxl by usixig satellite 
imagery, the ^)erBti«)ai sounder data are examiiMd as the soiree for estimating 
doud parametera as IniHit to the insdaticxi modetL Ih particular, seven of the 
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Dependent . ^ In^pendent Imagery HIRS. AVHRR 

Dataset I ^ Dataset * 
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FIGURE 2 BLOCK DIA6RAH OF STUDY PROCEDURE 












HIRS ctmmuQa w<rt intro^ioed Into a mult^iLB dlaorlainant aniQyais 
pr^raa daaicMd to asttaato aloud t^pa. Tha doud typa la usad to apacUy tha 
doud tramalttanoa. (Whan raady tat oparaUonal lyylioatlto, toa doud 
apadfUMitlon and tha pradpltahla watar darivad fhoa tha aoundar data aay ba 
taatad In tha i ^ agp aaa l on avpr e aaio n). Tha organization of tha analydaprooadraa la 
outllnad In Flgura 2. 

4^ Modd. 

4^1 Bado Excraadon 

On phjrdod. grounda It waa daairad to modal a Unaar azpraadon 
dnoa rndtlpla Unaar ray aa d on waa to ba uaau to Mine tha algorithm. At ai^ 
ghm Uma tha total diraot-pluB-difnjaa adir radiatlMi Ps Inoldant at tlM aurfhoa 
oan ba axpraeBail In tarma at aolar tlux Fq Inoldant on a horizontal aiarfioa at tha 
top of tha atmoaphara and an "aqULvalKit^ total flux tranamlttanca through ttM 
dkud-Araa or doudy atmoadMra. Ihdapendant imuiablaa are Introduoad through the 
tranamlttanoe, but preoautlon la taken to airald an axoaaa of variadaa that are not 
Indapandant of aaoh other. If T|^ r ap rea ants oloud-ftma atm<^3eharic tranamlttanoe 
and*^Q rap re aants dM oloudy atmo^haric transmlttanoa, than 

Fa*?o \ fo,]^(R/R)2coaZ (1) 

where 1 q la tha adar oonatant (1377 joulaa r'l), la tha mean aarth-aun 
dltotaiMa, and Z la the adar zenith an^. The Indloatad eoa Z aooounta only for tha 
horizontal aurfaoa at the top of tha atm opa h a r at dant paths Involving tha solar 
ztoith anda are impOidt in the transmlttanoa axpranwinnw. Eadi of tha 
transmittance terms oodd be fhotored ftrthar to separate, say, molecular 
scattering or absorption thorn aarosdaoatteringarabsorptlai. 

One method for retaining an exp raa d on for linear r^treadon with 
tarma that can ba asaodatad with phydcal procaaaaa la to conddar tha negative of 
tha natuml logarithm of Eq. (1) 
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r' ThiM, by d— Mm with tht log of ttos total tanM»altt«iiotg^o),«3totM*ixprwdon 
rtRilts In tnnifl of ttim lofKtthM of individual ataoaplMrio tnuMaittanoa taraa. 
To iKbaob that ttw tramaittanoa taraa oan ba r a pi^ aaa nt a d tv aj^onantiala, tha 
zl^it-taand dda c£ Eq. (2) rairaoaats ttia (^ittoaltii'ioknaaB. 

4^*2 YMiahHas and Paraaatariaattons 

(a) Nc3aou3ar Ataoapharat To aasoodata tha attanuation with tha 
rag^aaitnn ooaffloianta» tha doud-lhaa tranaaittanoa aiqr ba (^^praxtaatad by 

produot of twoaxponandalaa'^anda'^^^. Tha path langth Uiapoabulatadaa 
a praaara iwalad sdant path through tha total ataoqphario pradpitahla watv, aa 
aatiaatad fhoa tha airfhoa vapor preaaura a^ 

u a(bof bi iaXpa^Pr)M (3) 

where ^ oonstanta bo and are c.^aignad for an average r ap reee nt atlon of tha 
total preolpAtabila water (of. Smithl^ but oodld ba expanded to albw for aMoonal 
latitudinal variationB in the reHatlonehip to tha airfhoa daw point or» in thia oasa, 
the vapor preaaure. Tha ratio of atrfeoa preaaura to a refereroa preaatsvy aa a 
( aoaling fotbr, effaotivaly aooounta aa weE for large diangea due to atation 

elavation. Tha dant path magnification factor M r^ilacea tha aao Z ao aa to avoid 
problaffla with ttta plane^mrallal aaaum^ion at very low acan anglaa. Following 
Rogera, tha relaticxiahip ie axpreased aa 

M ■ k Qk2 - 1) ooa2 Z ♦ {] **i (4) 

with k 3 35. 

(b) Cloudai In order to maintain a aimpla initial dloud 

daaBlfLoation, each of tha aurtaoe obaervationa at oloudineaB ia placed Intc one of 
aevan total catagorleat CLEAR, HIGH, MID, FOG, LOW, MIXED and THICK. Tha 
oategorioal ranking ia in terma of an aaaumad inoraaatng atte n uati on (dacreaalng 
moud tranamlttanoa). Tha MID cmtagory actually Indludea h4di dloud alec, if they 
are p rea en b with tha mlddla olouds. Tha initial of FOG or 

OBSCURATION, LOW and MIXED were aubaequmOy altered on ttM baaia of 
tzmnamittancea inferred thorn tha SOLMET data and the initial reg reaa i oo mo^L 
The MIXED cat^ory initially iixsliuied FOG or LOW CLOUDS with MID and/or 
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S[Q8 doudi. Su b itq u i n Qy, aodof tte KZZ80 citoudi wnt nrtliwdfltrt aopordlm 
to tht doalMobliQrir dit«*miiiiag tiw ciloud otSdag (ftom UOom), with a adMqpMcit 
inortOM in tbo truvaittuiM for Vtm rtaainlag tOZBO oaUg ory . AJao, in tbo 
THICK o ata goty (tnniiidini nLnbooIzvtui «ad ouauloBlaburti tte oundooiabuB 
Mtofory was raolaadflad m low ftr aoatt«ad aB 0 UBt 8 » linoa tMr aflMa aay 
haft iMwa aMadL onr tha ohwrf to g atatlon. Obaounitiaa oondtttnni waia addad 
to tba THICK o a ta g or y, wbaaai LOW adaaqMOtay was a^vdad into ouaidua 
(CU)andflb«dtora C3T). 

To oonridar tba otloudy transaittanM in aora datail, tba 
inaolation aay ba d aac ri bad by 

P« • Fo [(1-C) \^C (5) 

whara C ia toa ciloud aaount and ia tba ciloud tranaatttanoa. For uaa with 
SOLHET data tba doud aaouxt ia d aa o t ib a d aa tba dsaarvad qpaque doud aaoia^ 
paua an ^ditlooal 35 p aroa nt of tba di f farano a betwaan tba total doud aaoiait and 
tba opaqua aaouat. Tha tranaaittanoa outdda tba doud di£f«a aoaawhab fhoa 
tba doud-Aaa banaaittanoa largaily as a result of aoattaring. Howavar, tbs 
dlfTerenoa can ba abanbad aa a faotor of <1* , parb^ps by lattiz« 'i* baooaa *. 
Than Eq. (5) baooaes 

Fa^o» % D-Cd-Oi*)] (6) 

Tba bradcat tara, raiaad to aoaa power (dose to unity) to aoooaodate tha 
approxLaatltxi and as a oonvanianoa for ragraaeton, daflnaa tb» cloudy 
treneafttanoa \and 

FyTo-Ta'Tc* aD-C('-'‘‘'*)]*2 (7) 

Tba aaauaad fora for is a Unaar fUnotion of tha coaina of tba adar zenith angLa 
with a tranaaittanoa factor, Y , to r a pr aaa nt tba tranaaittanoa for a partiodar 
oiloud typa at a particular solar zanltb angltt 

TCdo d<| ooa Z) (8) 
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kctaaOVt ttot cmnSL npr iBt ntatl o o of tht total tzaoiBlttonoa onr a Anita tlaa 
pmiod naad not naoaaaailly ba raotriotad to a liapia oloud "typa" with a gtTan 
aaount, aa iapBiad by Eq. (7). Ja fbot, tha cloudy tranaaittanoa tara ootdd ba a 
{Tottoot tora of aora than ona typa. Navathalaai, tar thda study it is assuaad 
toat tha oloud niaailfinafinn and tha tranaaittanoa ralattonship (an ba amiliad to a 
tlaa pariod of ona hoir. Tbaiafbrt, tha ooa Z that appaare in tha aodaiad 
STpraaslfinn is ooosidarad to ba rapi^raantad by tha awafs ooa Z for tha slvan hotr. 

(o) DaQy Araragasi Vhsn oonsldiring toa total ^dly inaolatloo, 
it baooaaa naoaasary to oocalda* tha avarafa daily ooa Z, whioh can bo axprassad 
in taras of too half-day boir angla H as 

oool^ «sln0ain6 * cob 0 cob 6 ain H/H (9) 

Hhara ^ is tha latituda and tha daOliintion. (Tha angidarvalocltyof tha aKth<>> « 
dh/db M 2 rad dajr'l). Tha aaxlaua solar zanith angla is glvan by tha dlffaranoa 
- 6), wharaas tha hour ani^ to aithar suorisa or sunaat. is givan by ooa H « -tan 0 
tan 6 . IT, during tha avaluatlon of Eq. (7) ovar a pariod of a day, tha nuaarator 
and danoainator of tha lafb-hand alda ara aua aad or intagratad sap^taly, than a 
raatrlotlon is iaplioitly lapoaed on tha avaraga d^ tranaaittanoa oorrasponding 
tothaoboarvad-ooaputadlaft-handaida. Ihpoitloular 


C A /V 

Fadt/V Fpdt* 



Ta^o^ldt (10) 

oSl 

iMbead of introduoing tha waicdi^ing fbotor indicatad in Eq. (10) aa^ of tha 
avaraged transaittanoas was apprcndaatad by asaignaant of a daily avaraga cos Z 
in pkca of tha 008 Z indioatad in thair definitions (Eqs. 4 and 8). 
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4.2.3 Regpaaalcm CoefUdente and Weighting Factora 

With the paraffletarlzations dtoecribed, the regression expreaal(xi 
for any hour becomes 


-to (F^Fq) = ao ♦ ai 0 - a2 to Tq 


( 11 ) 


whereas for the dally total 


^n (Fg^o) s -f S'] 0— a2to (12) 

where the caret Indicates that the dally average cos Z has be«i used Instead of the 
Instantaneous cos Z. to terms of the desired li»:ii}atiai Ita^, the houtiy or daily 
Insolation 


Fg s Fq exp (a© + a^ U - a 2 to tij)| (13) 

where the coefftoiaits have been determined by linear regression in terms of -to 
(Fg/Fo). 

When observations fbom several different hours (kj^) are to be 
combined in a sing^ estimate of the daily total insolaticxi) then additicxial 
r^res^n coefOcients are introduced in the analysis for each direct transmittance 

—to (Fg/F q) - a^ + a^ U + a2 In Tq (k^) + 

a3 to ^ (k2) + a4 to (k3) (14) 

wtnre chly a sin^ average U term is used for tt» day, but separate terms are 
maintained for the cloud transmlttances because it is not likely that they should all 
be wel^ted equally. Proper r^tlve wei^iting is handled by the re g ress i on 
coefadents. 

An alteniative procedue to that re p rese n ted by Eq. (14) for 
prespecified hours is to form a wei^ted average of separate daily InsolatLon 
estimates (uie for each observation time) based cxi infcx’mation obtained at 
differmt hours, characterized by polar-cxtiter overpasses. If Fg(ki) represents the 
estimated daily insolation based on data for the h(xr k-|. 
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where tile Wts rep re a e n t wet^riing flactore. More generally, if E ia uaed to define 
the estimated daily inaolation and tiie subacript k idoitifieB predictar hour, 

N 

E = E (Wit EicV j W|( (16) 

fc -1 

where aum mation is over ti» number N of daDy eatimatea baaed on individual 
howiy data. The wei^iting factor W)( wHI increase to a maximum far tbe hour 
with the daDy maximum cosine of the solar zenith angle, cos Z^x* 
r e pr ese n t s the time with the maximum corralaticm between the estimated E|( and 
tiM observed Fs» at least for tiM dependent data used in the r^^eaa l on analysis. 
Therefore, as a first approximation it is poaaLble to use the cotredatlcxi coefOciittits 
derived fnom ti» regression analysis as the wei^xting fbctor. However, the general 
form that was adopted is expresstd in terms of the solar zotiith angle for that hour, 
Z|(, as well as the maximum solar zenith angle 

W|(S A cosZmX'*’COsZ|( (17) 

where the constant A generally exceeds unity. 
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5U) RESULTS 

5.1 Obeerved Houriy-DailyaiacilatiMi CareiatiooB 


The SOLMET <tata fOeat with hotady insolation c^sservationa, are ideal 
for the (teterminatioa of the oocrdaticxi bstween an houriy total and the total daOy 
inao3athxi, the variabile of moat intereat. Such oorrsttattiona, apart finoa the ohvloua 
oontzibuticxi ariahig thorn the fhct that the houdy tc^al ia also pait of the daily 
toted, are especially informative about the relative aignifloarwe at obaarvation 
tlmea to be uaed in the daily eatiaation. Only fully auhUt hoirs are izvduded in the 
correlatians, which have been determined thorn data for 20 ditfemt stations 
dizdng each seasonal period. Biased (doud distributions are about the chly factor 
that coiild result in an aaym metry of the insolation distribution relative to solar 
noon. AEaf the results are summarized in Figire 3* Some di|ht skewness toward 
tlM afternoon hours is sqpparent. Othei*wise, the correlatioa curves are quite 
predictable, with the moat peaked diatributicxi for the shorter %finter days. 

It is aqq>ar«it that an observation near solar r»xn cairiea the most 
infcrmatiiMi about the daily totaL Obvioualy, a satellite obaervatic»i neer nocxi or 
eaidy afternoon is most uaefUL A lumber of heirs dose to rxxsi, say within three 
hours, also are bi^ily coerdated with the dally totaL The correlation ckxjps off 
near sunriae and sunset, everi if oedy full atailit hoirs are considered. 

Any ixMiividual staticxi skewness in distributiem is masked in Figire 3 
because each curve is based on data thorn 20 scattered stations. Seasonal 
variability in distribution spears largest in the aftemocxi. The ciorrelaticxis in 
Figis^ 3 m^ be interpreted as limiting conciltions for estimaticxi adnee they 
represent what would result if perfect hourly estimates were availabla (for any 
hour). To a first ai^proxLmatixu, the square of the correlatiiKB near sunrise and 
sunset also suggest the magnitude of the explained variance that might be 
associated with some skUL Thus, as wUl be seen later, re g re s sUxi results near 
midnight are assoedated with reduced explained variance that show littia value 
relative to the estimation of total dally insolaticn. 

5-2 Preliminary Commitatlcga 

5.2.1 Heuriy-Houriy Regressions 

Regression coefficients can be established betweoi hoiriy 
insdlatiun and hoiriy estimates based cn cdoud and water vapor slant path 
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FIGURE 3 DIURNAL VARIATION OF CORRELATION BETWEEN HOURLY 
AND DAILY TOTAL INSOLATION FOR TWENTY STATIONS 
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FIGURE 4 OBSERVED AND ESTIMATED HOURLY INSOLATION 
Group 1 - Bismarck, Dodge City, Madison 
Group 2 - Montgomery, Lake (Carles, Nashville 
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TABLE 3 

REGRESSION COEPnClENTS FOR PRELIMINARY 
HOURLY INSOLATION ESTIMATES 

(Standard emrs in paranttMses) 


GROUP 

INTERCEPT (a/0 

HATER VAPOR 
COEFFICIENTS (ai) 

CLOUD TRANSFER 
COEFFICIENTS (a^) 

1A 

.1765 (.0321) 

.0311 (.0051) 

1.2774 (.0338) 

1B 

.1059 (.0256) 

.0365 (.0032) 

.9771 (.0319) 

2A 

.1856 (.0283) 

.0272 (.0030) 

1.1240 (.0317) 

2B 

.1278 (.0242) 

.0278 (.0016) 

.8866 (.0251) 


1: Blsmwvk, Dodge City, Msuilacm 

As APR, MAY BsJUL, AUG 

2s Montgomery, Lake Chariee, NaAvUle 
As APR, MAY BsJUL, AUG 
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partmetarB. Figure 4 muatretes the initial aete of reaultal^ derived fhom two 
aeparate ffi'oiqps of 3 atatione eadi during aping and aummer nontha. (Gi^oup It 
BlaBarok, Dodge City, Madiaoni Oroig) 2s Montgomery, Lake Chaiiae, NaahvUIa). 
It ia apparent ftom theae reaulta (obtained without any fine tuning of initial 
panmetera) that there ia iitUa geogn^hloal dep«idmoe, but aeaacvial difreraxses 
ck> an^ear. Table 3 ahowa the ooneepo o ding regreaalon ooefOoienta and their 
atandard etrora. AH of the ooefOoienta are aigntfloant} magnttudea are 
re p reae n tatlve of thooe for other teata. CoefOeienta for the cloudy 
tranamittanoea, aa well aa intarcepta, are larger in the aping than the aum mo*. It 
waa anticipated that the cloudy atmosphere tranamittance coefOdenta would be 
doae to unity. 

The reaUlta in Flgire 4 Miow that the aummer eatimatea and 
diaervationa are aomewhat better correlated, with leaa biaa, than in apcing. The 
tdaa in the sping data ahowa overeetiffiatea for loweat magnitudea and 
undereetlmates for hi^er magnitudea. Part of the falUre to auooeaafUlly eatimate 
the extremes results merely Ihom the statistical nature of the coefOdenta. Some 
of the low-magnitude overestimatton also could reault from urxtarestlmateo of 
doud absorptton in oertair. douds, whereas in other drcumatances the doud 
tranamittance is underestifflated. Improper apectOcations of the i^iysLcal processes 
affecting the linear regreasLai modd. are likely sources of nonlinearity. On. the 
other hand, evoi with the proper regression expreasian in haiulllng the logarithm of 
the insolation, some distortion is Introduced in the average of the Inadation 
estimates that are baaed on the least>squares coeftldenta for the logarithm of 
insolation. 

5.2.2 Houriy-Dafly Estimates 

IT houriy estimates of Inadation are obtained throughout the day, 
tiiey may be aum med to provide estimates of the daily total insolation. Figire 5 
presents such sample results for cme group of stations shown in Figire 4, for both 
spring and aummer. Presumably, audi estifflates are auperlcr to those daily 
estimates that do not irxilude estimates for every hour. As can be sew in Figire 5, 
tlMre are far fewer paints («ily one per day) than in houriy-4iouriy comparisons. 
Fiarthermore, the averaging over the day results in a degradatioi of significance at 
the U coefficient, which is dependant cn the diimal variations in the dant path 
(the average daily zenith angle does not vary significantly within the data sample). 
As before, the ^ring estimates are somewhat more biased than the summer 
estimates. 
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FIGURE 5 DAILY INSOLATI(»i, OBSERVED AND ESTIMATED 
FROM HOURLY SUMMATIONS 
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FIGURE 6 DAILY INSOLATION, OBSERVED AND ESTIMATED FROM 
REGRESSION FOR COMBINATION OF PARAMETERS 
AT 0730 AND 1500 LOCAL TIME 
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If « 0 tiDfttes are not available for eadi hour, ao that the hourly 
eatifflates oan be aumned to obtain the dally total, the doud and aoiature 
obeervati(X)8 at any hour oan be used to eetimate the total daily insolation on the 
baaia of houdy-dally regreoaion (see Eq. 12). Or, if auitabiLe weighting coefflcLents 
w«« avaUabla, Individual houriy estimates could be wei^ited to provide a daily 
estifflate (see Eq. 16). Another means for obtaining dally eatimatas fhom hoixSy 
data is to use regreasLra ooefOeients that have bam designed to link parameters 
from two or more qpeoifLc hocrs direoOy to the daily insolation (see Eq. 14). For 
the first test of this latter procedure, water vapor and cloud parameters for 0730 
and 1500, approximating NOAA-6 and TIROS-N daytime passes, were used to 
specify the daily insolation. Regresaim coeftLclents for the doud transmittance 
terms are dominant but smaUsr than those associated with a sin^ hour. Results 
cf the preliminary oomputatims cocre^nding to the data in Figure 5 are 
illustrated in Figure 6. Deqpite a SU^ reduction in the correlatiai and the dight 
increase In residual differences the overall distilbutims of estimates and 
dsservations are similar to those for the hour3y summations. This result was 
encouraging, since the polar arbiters do not provide the lhaqumcy of observaticm or 
the optimum times of observation that are iiMduded in houdy summations. 
However, it does appear that the bias in ^xing is more pronounced than for the 
hourly sum mations, althou^ some overestimates of the daily inadation are now 
made at the hi^magnitude end of the range. 

5.2.3 Apparent Sensltlvitlas 

The preliminary computations indicated that the regression 
analyds approach was feaslbile, but that improvements in parameterization would 
be worthwhile. Regression coefficients for the ptsdpltable water term appear to 
be significant ordy when significant variations in the average cosins of the solar 
zenith angle occur. Surface dewpoint probably is a;i adequate estimatar of 
effective water vapor path. Clouds dominate the results (also rdated indirectly to 
surface dew point) and their parameterization is the moet lexical area for 
lmprx>vematt. 

Initial computatixx 0 for selected U.S. staticxv east of 105W 
longitude and between 30N and 47N latitude failed to show any pronounced 
latitudinal cr regional influence on the regression coefTlcienta. On the otherhand, 
seasonal variations appear to be significant, both in terms of the magnitudes the 
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rapi'tagdfEXi ooefflolantBi and the joint diabibution of estimatiraw and obeenratlonB. 
Hw distributions show aoae bias in tlMt astiaatasi e^)ecially with deputure fhom 
local nocxu 

The questiM of the dependency of results on individual stati«» 
was checked flirther by computing r^pesston ooeCTLeients separately for each of 20 
ataticms, aiKi ooeftloier^ also fbr aZL 20 stattons combined. Table 4 summarizes 
the moe; important regreaston ooefTiaient (for the doud transmittance term) for 
the most significant hour, 1200-1300 local time, appropriate to two 4-moibr 
periods, AprUi-May-October-November and June-July-August^ptember. The 
average coeffldant is largest In bie e^ring/falL period. Urdy two at the stations 
show anomalies of twenty percent or more in the same seme for both seasons. 
Fairbanks has low coefOoiaits airi Sterling shows anomaloudy hi^ coefficients. 

5.3 Modifications 

5.3*1 Regression Expression 

With the dataset expanded to 20 statlcxs and SOLMET data over a 
2l year period, checks van made on any improved correlation between estimatUxi 
and observation as a result of alteraticxis in the r^resslon e xp r es s i on. Some 
improvement was found by the simple introducttori of the doid amount C as a 
factor in the logarithmic expression (Eq. 12): 


•to (FjyTo) s So ai U a 2 C(Jn %) (18) 

Eq. (18) was used in aH subsequent r^resslon analyses. The improvement, while 
not explained on physical graui^ may have been related to the nonhomogeneity of 
doud ccxiditions (Le., d^)atture tpom an idealized single doud cover). 

The terms with carets in Eq. (18) indicate that a single average 
cosine of the solar z«iith angle is used to describe the path for the daily estimate. 
A comparison of results using ebs^, the daily average cosine, or cos Zmx» 
maximum cosine, showed very little difference in results. Consequently, most of 
tha final regressions that were performed used cos Z^x instead of cos Z. This 
change is justified in part by Eq. (10). 

The impact of introducing doud amount C in the last term of Eq. 
(18) is illustrated in Figure 7, which shows the distributinis for ertimates and 
obsttvations both without fr(a]) and with(7(b^the factor C. Ody the results Tram 
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TAsyLl 

SUMMARY OP DOMINANT DiSOLAYION-ESnMATION REGRESSION 
COEmCSNT COMPUTED FOR INDIVIDUAL STATIONS AT 



1200-1300 LOCAL TIME 



, CLOUD TRANSMITTANCE 
REGRESSION COEFnCIBNT 

STAHON 

SPRING/FALL 

SUMMER 

ChKi« 

1J)20 

0.746 

ColunUa 

1J)56 

0.804 

Miami 

0.992 

0.626 

BrownsrviHs 

0.800 

0.694 

RaLei^ 

1.174 

0.680 

Moabfpomary 

1J180 

0.586 

NaabvUlm 

1.119 

0.812 

Dodga City 

1J)13 

0.802 

Caribou 

0.873 

1.074 

Burlington 

1J)64 

1.017 

Madiaon 

1.035 

0.949 

Midland 

1.067 

0.963 

Biamarok 

0.804 

0.888 

Graat Falla 

0.865 

0.960 

Fairbanka 

0.499 

0.629 

StarUng 

IwJIO 

1.088 

TaBahusea 

0J106 

0.598 

Didianivoilia 

1.154 

0.809 

Pittaburg 

1.057 

0.818 

JOINT COEFFICIENT 1.000 

0.836 
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ttM 6«ti^lOR dst«Mt tor aurs a«r* art tfwwn. Tha oocraiatlon ocaCOoiaiA tor t2iia 
oaaa ia iaoraaaad IhOB 88 to 92 r«^^oant. 

5.3^ Cl oud Traaiatttaaot Panfatariaatlop 

For tiM OKllaot (prollalnary) ooaputaSiom idth a aaalL dataaat, 
ntgraarion ooaputatiom In a oo or dan oa with Eq. (12) uwd a ^van aat of oDoud 
oatafociaa^ aaoh with a ftvan Y fa* ooapuUne tha tnnaatttanoa fhoa Eq. (8) wi^ 
do and d^ aat at 0.6 and 0.7i re^paotivaiyt 

Y (0.6 ooa Z -t> 0.7) . 

Tba dxdoa of do and d^ raqidraa Vtmlt y la tha tnavaittanoa for ooa Z « 0.5* 
fiittially adopted valuaa of Y aralndioatadinthafltBtllatof Y^inTaUaS* 

Aftar tha pnillBlnary ooaputationB with tha initial oikud 
oatagoriaa and y % tha aotlBatae of inaolatian wara axaninad by oloud oatagcry. 
On tba baak of tha bias. If any, obaa?vad in aadi oatagory tha y ^ wara aitarad 
aooordin^. Rowavar tha o ona tawa do and di (Eq. 8) also wore aitarad ao aa to 
^dna Y aa tba doud tranaaittanoa for an ovarhaad aun (oos Z « 1.0). Thus, tha 
aaoond aat of tranaaittenoaa that appaar In Tatla 5 wOra aitarad both to 
ooapanaata for obaenrad Uiaas and to noraaUza to an ovarhaad aau Thla aaoond 
aat of tranaaittttwaa kloc^ with « Y (0.7 ooa Z •*> 0.3)) waa uaad for tha 
intoractilva data taata, ualng aataUite iaagary tor doud paraaatar apadflcattlnn. 

Coaputationa with tha second sat of doud tramaittanoas wara 
a^paratad into thraa oatagoriaoi CLEAR (do.1 oloud amount), CLOUDY (hO.9 doud 
amount), and tha REMAINDER. Moat of tha d^ fall in tha lattar oatagory. 
Ragraadon results for tha lattar oatagory wara dailar to thoaa dstainad 
pravioudy. Rapaadon andyaao for tha CLOUDY oatagory wara atrpridBg iti chat 
tha axplainad varlanoa wm vary low. As a result of tha ralatlvaly poor job of 
fitting tha CLOUDY oases, altarations were soirght for iaprovaaant of the 
paraaatarization. To ^tain a basis for a Bodifiad panaaatarization, tha 
ragranwlno ooaffiolaiits and tiM observatiom wara introduoad into tha ragraaslon 
expression. Solution for tha Y Is was addavad by invartiiy tha ragraadon 
expreasion. At ttm saaa tiaa, arfhea reports were jsed to divlda the LOW doud 
oatagory into alratifora and euauUfora. Mud) of the MIXED doik oa tagor y was 
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{daoaddjblierJntoaiMaftheLOW categarias or ^ MID catagocy, aooording to tiM 
observed eeOiiHt. TbersdeOned y ^ that ^^Msr as tte fixiil listing in Tsbls 5 were 
baaed only on a ata^ ooaputationaL abm> rallMr than an interactive procedure. 
They r tqr eae ob elfeotive empiric trananittancea, depenttacb cn the airfhce 
(^aaervaticsia, v^at were designed to yield laproved eatinatee of daily insolation 
ftoa ttm modi^suk^jted. 

Figure 8 iDustrates aeveral of tiw effective cilowl transaittances 
(total direct and diffuse) sa a fUncticn of tiMCQBine of the aoiiar zenith anc^ For 
coapariaon, aoae of the tiieoretical coaputations by Liou‘^5 are iiKduded in Figure 
8. mie large di8ca:vq)anc tetiiew transaittances for the cuaulus category results 
to aoae extent fifx>a ihet that the observations are avaflabls only fbr partial 
oloiid covers and include erhanoed aoattering fhoa doud waDa. 

Ih additioa to the diangea in doud oategoriea, aoae ot tte 
transaittances were increased by a aignifioant aaount. An objective ansQysb of 
the manner in whidi all of the doud catteries would be identified in the satdlite 
data was not completed. NevertbetLesa> the iapact of revised doud categories of 
the SOLMET data was to improve the insolation estimates. 

5.3*3 Final Adjustment 

DaasB\K!h as regresshxi ooeffldgrta were determined in a 
least squares procedire ai^slied to the logarithm of the insdation, it was . 
ocmddered likely that a statistical bias codd result in average urxierestimates of 
the insdation itself (when ai^dying the coeffidaits). A possible method to avoid 
such bias is to apply a ixxilinear least squares regreaslcxi analysis cUrecOy to the 
insdation fora of the equation (see Eq. 13)> in wMdi the coefficients all appear 
within the aiguaent of an expon«itial tera. In this approadi the attempt is to 
obtain coefficients that minimize the residual insolation standard deviatkxu 
Figure 9A shows the resulting distribution of estimated and observed daily 
insdaticxi based on nonlii^ar regressicxi with data for the best local hour, 1200- 
1 300. For comparison, Figure 9B shows the distributini for estimates derived after 
tte usual linear r^P'ession on the logarithm of the insolation. It appears that the 
noriUaear technique did not result in a superior estimate (slightly worse in this 
case). Thus it is likdy that any r^esidual estimation bias, if not fbom i^i^aLcally 
inadequate parameterization, mi^pit arise fbom the statistical assodatiicxi of data 
for any <X)e hour with that for the entire day. For example, if three separate hours 
are involved in separate estimatioiB of the daily insolation and aU Aow doudy 
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oonditlonB, tbm each eatimata wSI ba biaaed baoauaa aaoh regneiim wiU 

induda daar data fbr otbar bom rathar than conaidaring the ^laciflo Joint 
obaarvatlcxi at tha othar two aatlmatini houra. If treated aLmultanaouaily in a Joint 
r^pnaaaLMi, the fbct that aE three houa ware cloudy would indloata a much batter 
dmuaca for the oitira day b^ng cloudy and would raault in lower magnitude than 
fbr Vtm individual hourly treatmonta. SLmilaziy, the dngla aquatioci might give 
hi^MT magnitiKtoa fbr dear aky oo^dticxia ttian doaa the weighted aum of individual 
aetimationa. Theae pomthllltiiw auggeat that the wei^iting fhctora to be applied to 
iiKhvidual hourly eatimatea of daily inaoilatiCHi should take into account the known 
prevailing sky c<»ditlona at all (three) timea of daily eatlmatiini. 

If redduaL estimation b<as is aufflciaiUy large, it can be removed 
an approximate corectlni applied directly to the estimate for a givoi hour. This 
can be accomplidied in the following manner wtMre E|< is ttie original estimate and 
E|(' is the revised estimate: 


Ek'= 00^*01 Ek + C2Ek2 + 


(19) 


In effect, (xdy a sli^ linear rotation may be required. The 
coeffldents in (19) are determined by i^)ecifVing values of E|(' fK>m the observed Fs 
in the dependent datasets. Subsequaitly, the revised estimates of Fg, through Ei(', 
can be applied to the estimates that are obtained for the wei^ted combination of 
data fhom three different hours. Alternatively, Eq. (19) could be applied directly 
to a weij^ited estimate of the daily insolation. 


TABLE 5 

INITIAL AND REVISED CLOUD CATEGORIES 
AND TRANSMITTANCE FACTORS 



INITIAL* 


SECOND** 

FINAL** 

CATEGORY 

DESCRIPTOR 



DESCRIPTOR 


1 

CLEAR 

1.00 

.98 

CLEAR 

.99 

2 

HIGH 

.66 

.80 

HIGH 

.75 

3 

MID w/wo HIGH 

.50 

.63 

CU 

.70 

4 

FOG; OBSCURAnON 

.40 

.50 

MIXED 

.65 

5 

LOW 

.33 

.40 

MID w/wo HIGH 

.60 

6 

MIXED 

.22 

.35 

ST 

.48 

7 

THICK 

.11 

.22 

THICK 

.30 


• VrO.6 +0.7 
-0.7 +0.3 
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5.4 Matmottva Teats Fro« SatwIUto Photoemphs 

To oonduot a UmJtod "opamtlonal" tast of tha taohniqua for aallmating 
total daily Inaolatlon, haiM oo^ prints of reatrlotad aatalUta iaagery ware auindisd 
to Intan^ve analyst. Prints ware ganarslly available at raaoihrtlonp of about 
8km or lower, with three to four infimred images oonsidared aooeptable and with 
only a aii^ visible print allowed in the test. Use of the low raaoHution aing^ 
channal GOES infhared data was expected to ]deld estimates inferior to those 
determined flrom SOLMET data. NevaridMleas, the tests indloate the overall 
adequacy of the «q)proach. Of oounse, wh«i the model with the SOLMET data was 
formulated, it was not geared to the interactive analyst working with data at 
relatively low resolutim. Ihdependmt data for testing the daily insolation 
estimation technique were acquired for a spring period (15 April to 15 May 1979) 
and a summer period (August, 1979). Ten SOLMET stations were used, distributed 
in fai^ plains (Midland, Dodge City), the prairie (Bismark,Omaha, Columbia), the 
southeast (NaahviHe, Montgomery, Raleigh) and the Ohio Valley Qhdianapolis, 
Pitt^urg). For each station on each day, obeervations of cloud amount and ^rpe 
were made Ihsm GOES fhllrdlsk IR imagery at times approximating those of the 
two polar orbiting satellites (3 a.m., 7 a.m., 3 p.m., 7 p.m.). One eariy afternoon 
visible image was used as weQL In addition, dewpoints were recorded fl?om 
hemis(here surfhce maps for each observation time, for use in the preclpltablB 
water term. 

5.4.1 Spring Tests 

SOLMET insdatim data were oily available for 196 of the 300 
'^ticoMlaya^ for wMdi data were coUecsted. Three of the four obeervation times 
were used (eariy a.m., late p.m. and evening). Ctkch ^ Lae three observaticxB was 
implied to the particular predictive equation correqxxtdlng to the time the 
observation. This yielded three difTering predlctiixis* of the daily total insolation. 
The three values were combined Into a single estimate through weighting 


^Implicit in the development of each esti.mati(xi equation is the tacit aasumptlcxi 
that ccxuiitions obaerved at that hour wiU prevail throughout the day. A daily 
average of the cosine of the scQar zenith angpa is inciluded in the parameterization. 
Regression results appear to be relatively insnoitive to dianges in the efTective 
solar zenith angle included. 
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taohnlipMB (the sum of the products of oetiMetee and wdl^itiag flustore dl'vided by 
the aua of the weighting fhctor^. The wei|dttii« prooedire eei^Mins the valuos 
doeaeb to noon, the beet tine fhreetlBetes of daily totala. One to<dnique uead the 
ootreletlon oorffldenta, aaaoaieted with the iggrcgriete houriy equatione prochxsed 
in the develqpmmtal stage, aa the weights. A second tedmique aodifled moh of 
these weights by auUtipiUoation with the eddttionaL fiMtor 0 (cos Z/oos Z^xO 
where the cosine of the aolsr zenith an^ cos Z, refSre to the airarage for the 
pertioularhour whereas cos Zhx refers to the solar noon vaiiM. 

As an additional test, use was made of a sii^ predictive 
equation with ragresdon ooefflciMits baaed on the incorporation of all three 
dservatlma at the three different times in one regresBdon esgresalan. Only one 
estimate of the daily total insolation is obtained, and no weitfiting is involved in 
this technique. 

5.4.2 Sum mer Teebs 

The summer tests were oonduoted in the same manner as the 
spring tests, and with shout the same ratio of actual data matohea to the total 
numbers of estimates (198 to 310 atatiMi days for summer). Ih additina to the 
1 thrae weighting methods (combining daily astimates firom three aeperate 

obeervation times, and the ainglfi predictive equaticxv), a third weighting method 
was introduced for the sum mer tests. This last weighting method was intended to 
be fl?ee of the ^)eciflc ragreaaion results and thus did not involve the cocrelation 
coefficients. It did izKilude the cosine of the maximum sdaT zenith angle (a 
functioa of latitucte and declination} and the cosine of the solar zadth angle for the 
hoir of estlfflation. Therefore, the noo-cotrelative weighting method employed the 
weighting factor 

jj.1 cosZmx ‘•'COszJ 

In this expreesicxi the cosine of the solar zonith angle could be allowed to become 
negative, since the coefficient of cos Zmx would prevent the total wei^iting 
factor ihom becoming negative. 

Dudng the interactive analysis for the summer season it was 
posBlhle to introduce significant errors in tranaitioa fhom a satellite image to a 
representative parameterization of dcudineas at a local point. Most of the thick 
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doudizMM was of tha ooivaotlva ^pe with ocMMUkrtbla variation ovar abort 
diatanoea. This may bava oontzibutad to aoma dagnadatian in quality of tha 
aun fflar raaulta as oompared to those of spring. 

5.4.3 Raauilta From Meraotlva Tasta 

Table 6 auafflarizea the qpring and aummor taA results obtained 
with tha different weighting mathoda. la both taata tha t fa r aa - e q uatioii (indlvidittl 
hourly res^eaalona) modala with weighting of each dally aatifflata peiformad batter 
than tha single aquation model (single ra gi ^eee t on for oombdnation of houra), at least 
in terms of corrailation and tto atandard deviatlcm of the reeiduaL The beat result 
in both testa sterna lh>m ttM modified cocrelation wei^iting procedure, but a teat 
with the nonp«orrelatlva weighting stK>wed about the aame results 

Table 6 

RESULTS FROM INTERACTIVE TEST WITH SATELLITE IMAGERY 


WEIGHTING METHOD 

TEST 1 (CotreilaticMi weighting) 

TEST 2 (Modified cocralaticMi weighting) 
TEST 3 (Single pridictlve equathxi) 

TEST 4 (N<M>carrelatlva wei^iting) 


CORRELATION 

RATIOi ReaidualRMS 

EstLmated-O baerved 

To Observed Mean 

Spring 

Summer 

Spring 

Sum mar 

0.903 

0.819 

0.191 

0.184 

0.904 

0.836 

0.185 

0.178 

0.888 

0.794 

0.239 

0.184 


0.831 

■ - 

0.179 


m general, it was anticipated that olearMiay insoilatiai would be 
underestimated aomewhat, while 61oudy-day inaolaticai would be overestimated. 
For the interactive teats, the second or intermediate set of doud categories airi 
transmittance factors was used (see Table 5). Some difficulty was expected with 
an Interactive "estimate” of doud amounti sudi estimates usually invdve a 
regk»ial averaging because of the im practicality of pinpointing a usefUllocal cloud 
amount lh>m a fbH disk-image. With the imagery available, dx>ud amount was 
qpadfLed (xlly to the nearest quarter of the aky. Nevertheless, even thou^ 
c^mum data were not used and not all of the recmt model improvements were 
included, the aum mary In Table 6 indlcatea that the interactive testa were quite 
atnceeafd. 
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FIGURE 12 OBSERVED DAILY INSOLATION AND ESTIMATES BASED ON SINGLE PREDICTIVE 
EXPRESSION INCORPORATING DATA FROM THREE DAILY SATELLITE IMAGES. AUGUST 1980 



Flgura 10 iSuBtritas results for tiw spring test airi Flgire 1 1 
abows tbs auaaer resUltB. AUhai;gb tbs spring nesults abow audi 1mm bias than 
thoss of thsauBBertsst, andabow Ugberooire]atloo,tbsrmtioof tbssRvrtotbs 
^saenrsd BMn is aobaOy a bit sBaUsr for tbs BiBBsr test. With ths Mgber 
auBBSitias Bean this is so dsqdts Vtm Bore obvious bias in tbs suBBer resUlta 
(aoBswhat sxaggerstsd hir tbs soaLss). Figurs 12 abows another estlBsts for ths 
auBBsr tsst data, this tias using tbs tbres hour ooabinsd ro^asaion squation 
Ostn^ prsdiotivs squation). This result abows Isss Wessfi sstlaatss On ths 
sxtrsBss) than tboss in Figure 1 1, altbou^ \3am oorrelation is lower and ths R MS a 
littls larger (aes Tabls 6). Not abown is tbs fact that ths averags dsvlation 
(absoluts) of estlaats fhoa observation was actually SBallar for ths ain^ 
prsdlotivs squation. Future dsvel^aent ateuld rssxaBins the posaibility of 
•djuatii^ ths aingls squation approach, with Bodiflcation, as an altemativs to 
iaprovixig ths weitfitiiig fhotOTB for individual 
5.5 Final ReareariCTi Modtol Rswflts 

SOLMET date were apiUt into tte«s groups for final analyaiat 
DsosaberWanuary-Fsbruary-MardiCwinter), Apidl*May^}ot^)et«-NovoBber(aprii^ 
fan), and Jui^July^Auguat-SsptSBberOaunBer}* Ths qpring-AU ooabination 
aubsequoitly was aeparatsd into ^ring and fall subgroups, each about half as large 
as tbs winter and auamer ^cups. Attasuatiai reffmdm coefficients, along with 
tbsir standard errors, are listed for eadi hour the day in Tablae 7»4. Tlw 
nuBber of observations Involved in earit hoiB^yi^ally analysis showed a alnlBUB 
after aidn4d>>t and a aaxlaua in aid-aftsmocn for all s e as on s . & teras of tlM 
esUaated insolation itself, after exponentiation with the derived cosfflcittits, the 
correlation ooeffloients between estiaatad and ^jservsd insolation are dsterainsd, 
as are ths standard dsviatlm of ths reeidiad dlfrerencea. These teras are included 
in ttw tables. 

Seasonal Granges in ths estiaation observation statiatlcs are evident. 
Correlation coeftldants, disregarding the nighttlae values, ar« audi thr aaae in 
qxing and suaaer, h4dMSt in fbll, and alaost as hi|^ in winter. Ths standard 
deviation of reriduaL differences, or standard residual, abows a general inverse 
relatiradiip to the correlation. It was saallest (bring ths fSU and winter. During 
tbs iaportant qring and aua aer seasms, the ratio of the Aandard residual to the 
asan ^»ervsd insolation bscoaes saaDer in suaaer with the higher observed 
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TABLE 7« 

SPRINa 

SUMMARY OF DAILY XNSOLAISON COEFUCSNI^ DETERMINED FROM 
RBQRBSaON ANALYSIS FOR AVERAGE BACEOROUND ONTERCEPT), 
WATER VAPOR SLANT PATH, AND CLOUD ATTENUATION PARAMETERS 

Stindird Mttra trt in ptrsntiMsts 


HOUR 

NUMBER 
OF OBS 

ZN^CEPT 

(»l) 

PATH 

(A2) 

CLOUD 

CORR 

01 

641 

a72l4 (0U)57) 

-0.0205 (0X19) 

0.5505 (0.054) 

0.46 

02 

641 

0.6706 (0U>56) 

-0X185 (0X19) 

0X116(0.055) 

0.49 

03 

1141 

0.6064 <0U)41) 

0X089 (0.014) 

0.6991 (0.042) 

a53 

04 

1187 

a5872 (0U)39) 

0X128 (0X14) 

0.6967 (0.040) 

0.55 

05 

1187 

0.5409 (0X39) 

0.0247 (0.014) 

0.7417 (0.041) 

0.56 

06 

1356 

0.4771 (0X34) 

0X335 (0.012) 

0X127 (0.037) 

0.62 

07 

1502 

a4508 (0X31) 

0X278(0X11) 

0.9090 (0.034) 

0.68 

08 

1464 

0.4118(0X30) 

0X249 (0X10) 

1.0642 (0.034) 

C.74 

09 

1487 

0.3654 (0X27) 

0X323 (0.009) 

1.1678 (0X3D 

0.79 

10 

1521 

0.3390 (0.025) 

0X345 (0X09) 

1X456 (0.029) 

0X3 

C': 

1504 

0X965 (0X23) 

0X377 (0X08) 

1.3682 (0.028) 

0X6 

12 

1549 

0X666 (0X22) 

0.0436 (0X07) 

1.4058 (0.027) 

0X8 

13 

1585 

0X427 (0.021) 

0X497 (0X07) 

1.4344 (0.026) 

0X9 

14 

1575 

0X047 (0X20) 

0X603 (0X07) 

1.5110(0.026) 

0X9 

15 

1591 

0X176 (0X21) 

0X671 (0X08) 

1.4454 (0.027) 

0X8 

16 

1604 

0X442 (0X22) 

0X667 (0.008) 

1.3848 (0.027) 

0X6 

17 

1556 

0X584 (0X24) 

0X660 (0X08) 

1.3701 (0.030) 

0X3 

18 

1541 

0.3205 (0X25) 

0X525 (0.009) 

1X659 (0X32) 

0.78 

19 

1546 

0.3853 (0X27) 

0X389 (0.010) 

1.1641 (0.033) 

0.74 

20 

1417 

0.4269 (0X30) 

0X298(0X11) 

1.1182(0.036) 

0.70 

21 

1332 

0.4744 (0X33) 

0X230 (0X12) 

1.0388 (0.038) 

0.67 

22 

1318 

0.5208 (0.034) 

0X167 (0.012) 

0.9693 (0X38) 

0.64 

23 

1179 

0.4394 (0X37) 

0.0670 (0X10) 

0.9018 (0X45) 

0.59 

24 

746 

0.5956 (0X44) 

-0X037 (0X16) 

0.7598 (0.049) 

0.58 


HOUR s Looiltlmtatandaf hor 

CORR s Con^slatlon b«tw«m aatimatad and obMTvtd daily inscilatlon 
STD RES « Standard daviation of raodduals (Kj b~2) 


STD RES 

6620 

6m 

635M 

6268 

6182 

5912 

5537 

5074 

4583 

4170 

3784 

354^4 

3341 

3356 

3595 

3871 

4239 

4695 

5048 

5414 

5691 

5891 

6233 

5998 
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SOMHIR 


s 


t 

SUMNARX or DAILY DISOLA1ION COBFFICIBNTS DETERMSIBO FROM 
REGRESSION ANALYSIS FOR AVERAGE BACEQROUND CENTER CEPT), 
WATER VAPOR SLANT PATH, AND CLOOD ATTENUATION PARAMETERS 

Stindard mm art in pTtntluBM 


HOUR 

NUMBER 
OF OBS 

IN^CEPT 

(ai) 

PATH 

(M9) 

CLOUD 

CORR 

STD RES 

01 

1325 

0.4217 (0.043) 

0X556(0X11) 

0.4996 (0X32) 

0.59 

5610 

02 

1325 

0.4136 (0U)43) 

0X557(0X11) 

0.5279 (0.032) 

0X0 

5563 

03 

2197 

0.5189 (0X31) 

0.0329 (0X08) 

0.5885 (0X26) 

0X0 

5664 

’04 

2283 

a5035 (0X30) 

0X345 (0X08) 

0.6019 (0X25) 

0X2 

5535 

05 

2289 

0.4669 (0X29) 

0X403 (0X08) 

0.6617 (0X25) 

0.64 

5372 

06 

2752 

0.4114(0X24) 

0X421 (0X06) 

0.7400 (0X22) 

0.70 

4970 

07 

3045 

0.4270 (0X23) 

0X308 (0X06) 

0.7751 (0X20) 

0.73 

4796 

08 

2961 

0.3918(0X22) 

0X317 (0X06) 

0.9174 (0X20) 

0.76 

4509 

09 

3028 

0.3676 (0X20) 

0X311 (0X05) 

1.0081 (0X19) 

0.8C 

4156 


3077 

0.3533 (0X18) 

0.0285 (0X05) 

1.1042 (0.018) 

0.83 

3803 

^ , 1 

3049 

0.3299(0X17) 

0X284 (0X04) 

1X350 (0X18) 

0.87 

3436 

12 

3203 

0.3074 (0X15) 

0.0311 (0X04) 

1.3056 (0X16) 

0.89 

3160 

13 

3254 

0.3159 (0X15) 

0.0312 (0X0^0 

1X673 (0.016) 

0X9 

3164 

14 

3225 

0X856 (0X15) 

0X376 (0X04) 

1.3304 (0X17) 

0X9 

3094 

15 

3274 

0X744 (0X15) 

0.0436 (0X04) 

1.3039 (0X17) 

0X6 

3247 

16 

3283 

0X975 (0X16) 

0X422 (0X04) 

1X235 (0.018) 

0X6 

3509 

17 

3173 

0X994 (0X17) 

0.0454 (0X05) 

1.1911 (0X19) 

0.63 

3802 

•18 

3156 

0.3318 (0X19) 

0.0457 (0X05) 

1X578 (0X21) 

0.78 

4296 

19 

3152 

0.3832 (0X21) 

0.0395 (0X06) 

0.9447 (0.022) 

0.74 

4704 

.20 

2823 

0.4103 (0X24) 

0X388 (0X06) 

0X952 (0X25) 

0.69 

5072 

21 

2706 

0.4239 (0X26) 

0X426 (0X07) 

0.7785 (0.025) 

0.66 

5297 

22 

2692 

0.4661 (0X27) 

0X372 (0.007) 

0.7027 (0.025) 

0.64 

5483 

23 

2331 

0.3918 (0X29) 

0.0665 (0X07) 

0.6531 (0X29) 

0.60 

5718 

24 

1589 

0.4637 (0X33) 

0X432 (0X09) 

0.5447 (0X30) 

0.61 

5648 


HOUR B Local tl 0 « stand of hour 

CORR s CoRvlatloa hetwocn ofltlBatod and djaarvod daOy izBolatlon 
STD RES s Standard dovlatlon of reddualaOO 
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TABLE 7o 
FALL 

SUMMARY OF DAILY INSOLATION COEFFICIENTS DETERMINED FROM 
REGRESSION ANALYSIS FOR AVERAGE BACKGROUND ONTERCEPT), 
WATER VAPOR SLANT PATH, AND CLOUD ATTENUAHON PARAMETERS 

Standand enrora are jn perentiieaee 


OUR 

NUMBER 
OF OBS 

01 

603 

02 

603 

03 

999 

04 

1034 

05 

1003 

06 

1111 

07 

1311 

08 

1292 

09 

1395 


1417 

, i 

1385 

12 

1447 

13 

1471 

14 

1448 

15 

1474 

16 

1493 

17 

1420 

18 

1378 

19 

1372 

20 

1221 

21 

1165 

22 

1169 

23 

1078 

24 

651 


(ao) 

INTERCEPT 

0.6111 (0.071) 
0.5886 (0.070) 
0.6279 (0U)55) 
0.5751 (0J)52) 
0.5550 (0.052) 
0.5355 (0.045) 
0.4960 (0.039) 
0.4463 (0.039) 
0.3868 (0.032) 
0.3586 (0.031) 
0.3210 (0.029) 
0.2882 (0.026) 
0.2858 (0.025) 
0.2787 (0.025) 
0.3010 (0.026) 
0.3395 (0.027) 
0.3505 (0.030) 
0.3624 (0.034) 
0.3989 (0.036) 
0.4403 (0.042) 
0.4829 (0.046) 
0.5656 (0.047) 
0.5366 (0.055) 
0.5807 (0.063) 


(ai) 

PATH 

0.0583 (0.019) 
0.0602 (0.019) 
0.0408 (0.014) 
0.0477 (0.013) 
0.0499 (0U)13) 
0U)462 (0J312) 
0.0435(0.011) 
0.0423 (0.010) 
0.0394 (0.009) 
0.0367 (0J308) 
0.0394 (0.008) 
0.0442 (0.007) 
0.0449 (0.007) 
0.0473 (0.007) 
0.0476 (0.007) 
0.0458 (0.008) 
0.0500 (0.008) 
0.0562 (0.009) 
0.0540 (0.007) 
0.0494 (0.01 1) 
0.0438(0.011) 
0.0303 (0.012) 
0.0721 (0.011) 
0.0482 (0.017) 


(a2) 

CLOUD 


0.4744 

(0.050) 

0.5154 

(0.051) 

0.5958 

(0.041) 

0.6383 

(0.038) 

0.6697 

(0.038) 

0.6952 

(0.034) 

0.7789 

(0.030) 

0.8522 (0.030) 

0.9590 

(0.027) 

1.0185 

(0.025) 

1.1035 

(0.024) 

1.1678 

(0.022) 

1.1818 

(0.021) 

1.1961 

(0.022) 

1.1537 

(0.023) 

1.0989 

(0.023) 

1.0811 

(0.026) 

1.0184 

(0.028) 

0.9543 

(0.028) 

0.9042 

(0.032) 

0.8298 

(0.034) 

0.7652 

(0.035) 

0.6410 

(0.041) 

0.6052 

(0.048) 


CORR 

STD RES 

0.62 

4310 

0.63 

4271 

0.65 

4082 

0.69 

3913 

0.70 

3786 

0.75 

3713 

0.79 

3425 

0.81 

3177 

0.85 

2821 

0.89 

2461 

0.91 

2155 

0.93 

1927 

0.94 

1869 

0.93 

1941 

0.92 

2133 

0.90 

2374 

0.87 

2635 

0.85 

2917 

0.83 

3105 

0.79 

3293 

0.77 

3421 

0.75 

3573 

0.69 

3913 

0.68 

3926 


HOUR s Local time at erKl of hour 

CORK s Correlatiai betwe«i estimated and observed daily insolation 
STD RES s Standard deviatim of residuals (Kj m*^ 
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TABLE 7d 
WINTER 


SUMMARY OF DAILY INSOL AHON COEFnCIENTS DETERMINED FROM 
REGRESSION ANALYSIS FOR AVERAGE BACKGROUND (INTERCEPT), 
WATER VAPOR SLANT PATH, AND CLOUD ATTENUATION PARAMETERS 

Standanl €nron «« in parentheses 


NUMBER 

(*o) 

(ai) 

(S2) 


OF OBS 

INTERCEPT 

PATH 

CLOUD 

CORR 

1658 

0.5878 (0.035) 

0.1362 (0.013) 

0.3258 (0.032) 

0.56 

1658 

0.5705 (0.035) 

0.1349 (0.013) 

0.3546 (0.033) 

0.57 

2800 

0.5968 (0J)27) 

0.0906 (0.009) 

0.4197 (0.025) 

0.62 

2880 

0.5688 (0.026) 

0.0887 (0.009) 

0.4592 (0.024) 

0.65 

2812 

0.5259 (0U127) 

0.0929 (0.009) 

0.5118(0.024) 

0.67 

2993 

0.4865 (0.024) 

0.0988 (0.009) 

0.5403 (0.022) 

0.70 

3251 

0.4407 (0J322) 

0.0984 (0.008) 

0.5903 (0.020) 

0.74 

3227 

0.3873 (0.021) 

0.0938 (0.008) 

0.6832 (0.020) 

0.77 

3490 

0.3315 (0.018) 

0.0894 (0.007) 

0.7745 (0.018) 

0.82 

3552 

0.3170 (0.017) 

0.0851 (0U)O6) 

0.8169 (0.017) 

0.85 

3461 

0.3017 (0.016) 

0.0781 (0.006) 

0.8978 (0.017) 

0.88 

3520 

0.2748 (0.015) 

0.0815 (0.006) 

0.9527 (0.016) 

0.90 

3571 

0.2586 (0.014) 

0.0814 (0.005) 

0.9963 (0.015) 

0.91 

3512 

0.2513 (0.015) 

0.0850 (0.006) 

1.0184 (0.015) 

0.92 

3582 

0.2563 (0.015) 

0.0923 (0.006) 

0.9907 (0.015) 

0.91 

3623 

0.2795 (0.015) 

0.0959 (0.006) 

0.9340 (0.016) 

0.89 

3492 

0.2901 (0.017) 

0.0970 (0.006) 

0.9078 (0.017) 

0.87 

3470 

0.3180 (0.018) 

0.1007 (0.007) 

0.8362 (0.018) 

0.84 

3476 

0.3636 (0.018) 

0.0946 (0.007) 

0.7803 (0.018) 

0.82 

3124 

0.3928 (0.021) 

0.0906 (0.008) 

0.7509 (0.020) 

0.78 

2982 

0.4366 (0.023) 

0.0861 (0.008) 

0.6923 (0.021) 

0.75 

2996 

0.4847 (0.023) 

0.0752 (0.008) 

0.6492 (0.021) 

0.73 

2825 

0.4983 (0.025) 

0.0688 (0.007) 

0.6368 (0.023) 

0.71 

1882 

0.5209 (0.028) 

0.0688(0.011) 

0.5385 (0.028) 

0.69 


HOUR s Local time at end of hour 

C 0 R R s Corr^tion between estimated and observed daily inscQation 
STD RES s Standard deviation of residuals (Kj m~2) 


STD RES 

M673 

4645 

4460 

4398 

4^1 

4123 

3913 

3662 

3290 

3035 

2729 

2470 

2314 

2280 

2410 

2646 

2902 

3221 

3435 

3619 

3763 

3893 

3985 

4078 
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The ataodard etrora of ttie reff^aaakm ooefOoienta are included 
parentheticaliy afbar eadi ooeffSoient in tte taMae. Only the few of the nltfittiDe 
water vi^xr path ooeffidnita durii^ aping (especially the n^ative valuea) are 
inatgniflcant. Otherwiae, the daytime water va^xr path ooefflaienta are amaUeat in 
aummer and largeab in winter. The moat important coefficient, for the cloud 
tranamittance term, ia largeat in spring and decQinaa with aeaaon to amaUeat valuea 
in winter. Ihtercepta taxi to be amaUeat in spring and aum mer and larger in faE 
asoA winter. 

AH co^ndenta vary diurnaUy. Near midnicpit, the explained vailancea 
and correlaticxi cofflcienta nada low levela, in many inatancea droj^dng below the 
variance aaaociated simply with the changing aolar zodth an^. Below auch a 
threehcld it ia COTsUiered that the re gr e a a l on rdatiauh^a for that hour coitain 
UtOe uaehsl information for eatlmatlcxi of the totel daily inaolatkxi. 

Figurea 1^-16 iUuatrate distributlcxis of estimated and obaerved daily 
insolation, baaed <xi tlM r^p'eaadon results that are available for the hour 1200-1300 
during each aeascxial period. Theae iUuSrationa give a good view of the scatter and 
residual bias in the data diatributiixB f<x* ttda optimal hour. Ih application, the 
daily eeHmate would be baaed on a combination of eatimatea fhom three differait 
timea. A distribution of such results would resemble the distributions illustrated. 

5.6 lndep«ident Data Teats 

5.6.1 IMepenctoibSOLMET Data 

Two months of SOLMET data, for an April and a May, were 
withheld tma the final regreseion analysis (depeixi«xt data) so that tl^ could be 
used as a aoirce of indspo^dent data for checking resulta. Two testa were made 
with these data, by using information taken fh>m three different hours (simulating 
data coverage fh>m two pclarKrbiting satellites). The <xily differoice betwe«i 
tests was a one-^x>ur shift in the time of the central data input, fh>m the 1300-1400 
local hour to 1400-1500. Ui each case the results were weighted in accordance with 
Eqs.(l6) and (17). For ttMse tests the latest or "final'* version of the cloud 
categories and transmittance factors (see Table 5) was used. For orw of the tests, 
an adjustment to remove bias, in accordance with two terms of Eq. (19), was 
attempted with factors based on the empirical dependent data for^^xing. 
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FIGURE 13 OBSERVED DAILY INSOLATION AND ESTIMATES BASED ON REGRESSION 

FOR 1200-1300 LST, SPRING 


LEGENOt « - t MS. • . 2 oSS. EfC 


ORIGINAL PAGE IS 
OF POOR QUALITY 


• n 

* • 9 ? ^ < 




o 0 8 < < • ^ 

CisSHK-'s -S 

• a'^BSssssSSss- 


«0oauik*»oe88 8 « 

<<»uuvii(8«SSt««a>S3««2* 

s-ladpsfggSIsSS*-- 

S S * 5 O IL S (ilHMli S * < • 

*: s,Ss 2 |siif*Hi :„** 

g9SUMHiuSzhS*S<SB « 


00 

ON>^ in 
OOvO iH 

• l-l . 

o<no 

I I N 


aasBssisSiSrsjis < 

**r5:sig839s- " 

: = l-fiiiiiiii!!. « . 

“ 5;i||82li|ls;i:5 ** , , 

« «« 35 *SsaB£a 2 <uw 3 < 

3 *:;i:SssHtaHS:*jj -*>* 

* S*2*2S299;i|S232S - 

- <9^ 239225592|. 2* 

a <• 22«8S Sg <5:^ 

- 5-< irs^Piipif 5^- 

^•22*92 25S88s^2,_ 


: 2 s 53 - * 

*<2 2S," "« 

^ 2 *-.S“ •' 9 a *'252 •* *2 
* «- •‘*8"«82 32*3, * * 

** ** 2*«s5i«l2l2*S2S3-92 
*a 9 32 a 22 « 2 ,* ,-2 


Z* -• 2 - 8 -*'- 5 '‘<*‘ < 

o**a* * 22 "* • 

• *a- 2 ** 2 *r: 2 -*S*- 

* •« S“* 2222 « 2«25 

* <a ,* «*5 * 

* , “!2 

< * « a a 


2 e 

o 


“ z 
o 

§ H 

« I 

f 5 

m* 

O 

I i 
5 2 

M 

a 

f CO 

I M 




1 S 

2 « 


I ^ 

9 8 


(j_® r^) NoiiviosNi aaAHasao 


50 


LCSCNO A - I ns. m m t o«, ttc. 


ORICSNau i ao. 1 - 
OF POOR QUALITY 


o\ 

f-4 

• 00 • 


I I I 


sag 


3 u < 

II** *• 

* 

9 9 a 4 


i9i** *•» 
««« 


/•ii:;’ nil; :« 

* 5*5- * 

• * S585S5 u 


i?s <• 


* S585S5 u , 

U: 

5*-5“I5s* •! 


•«*ogw5«<<e « «< « 

« **'‘S‘'w.<if*a.,« • 

“< «<®.3IS88 u*S“**« 

•JS-83S*S3r.S?u'‘ 

< j.SSfSSS^SIS 

-A'^sasossea 

^*5 *38838 

^ 4 <UO 0 a < 
^ <* 0 9Ca < 


T • 

I 2 


(j.® PI) NoiiviosNi aaAHasao 


r^T-^- 


estimated insolation (kj m'^) 

FIGURE 15 OBSERVED DAILY INSOLATION AND ESTIMATES BASED ON REGRESSION 

FOR 1200-1300 LST. FALL 
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FIGURE 16 OBSERVED DAILY INSOLATION AND ESTIMATES BASED ON REGRESSION 

FOR 1200-1300 LST, WINTER 


Reeults for the tests with Indepmdmt SOLMET data are shown In 
Figure 17a, b, and o. ]h the statlstloal sense It must be ocxiduded that the tests 
were oertalidy satiSTaotory, even thou^ satelUte Inputs were not actually used 
here. The shift In the midafternoon data souroe by one hour had a all^ not 
signifloant impact on the results. The adjustment to remove the 
re 0 ?eeBlon bias improved the RMS 1^8 percent aikl visibly improved the symmetry 
of the dlstiibution. The bias ooRsctim inre was linear, for the three constituent 
Ixan. A im-llnear or step^tinction correction mi^ improve the results further. 
5.6.2 ]hdep«ident Dlgtbal Jaftmred Sateilllte Data 

The aim of using digital miiltiqpeotral data fhom an operational 
pQlar-orbitiiig satellite is an automated technique for Insolaticwi estimatim 
involving doud attmuation Information. The Initial step of q)ecifVing doud types, 
for the parameterization of doud transmlttaiK^s was designed and tested as 
(Hacuased in the Appendix. 

Of most importance for an automated operaticxial procedure is 
the availability of an objactive doud type apedfk^tlon. Data from seven channels 
of the operational HIRS sounder were selected for analysis toward this 
deveLe^ment. It was assumed that high reeolutlon AVHRR data mi^ be used with 
a threehdding method to define a corre^)onding doud amant. Ultimately, both 
doud type aixi amount for a given resolution (say 40 km) could be passed along to 
the interactive analyst, who could assimilate this information with whatever 
imagery was available for making final estimates. 


53 


3n00a 


OWOWM. 

OF POOR QUALITY 


O 

r 



«f « «r < 


« 

e 




o 

e 


o 

u 

e 


9 

9 

O 


9 

9 

? 


9 

S 


s 


o 


•• ♦ • 

c 


o 

o 


s 

St 


^3-® r®) NoinosNi aaAHasao 

54 


I 


iS 

I ^ 


! CM 

( c I 
I c 

tv: 


S 

1 


- « 

e H 

,'S 

- o 

CO 

!? *-• 

j- o 

l.l 

- * 

I H 

I a W 

I 8 U 


IS 

♦ o 

I c 

I - 


\ t 

\i 


. ♦ — — — ♦ f 

p 


s 

1 

CO 

& ■ 
W M 

H0tf 

CO Oi 
U CO 

o • 

II 

iS 

H Ho 
W o 

o u 

2 CO - 
<«o 
O O 

Si-3 
nri - 
•<So 

»J oo 

O 00 

§e° 

M W-«^ 
O cd 
>• 2'-' 
I-} M 
H 0U 

SS 

QM 



s U 

8 

06 U 

bd S 

O 

• 

CO 2 


MH 


D 

U 



« • t m$. * ■ 2 nnSf ffc 



-9 


i $ 




9 


j 

H 


i 

i 


9 

9 


g 

m 


C 

<9 


g 


9 



(j.o m) NouviosNi asAHasao 


55 


FIGURE 17 OBSERVED DAILY INSOLATION AND WEIGHTED ESTIMATES BASED 
ON THREE INDEPENDENT DAILY OBSERVATIONS, SPRING 1979 

(b) 0800,1500,1900 
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6.0 CONCLOgONS AND RRCOMMENDATIOMS 

RMuitta (Maintd flron amlyait of SOLMET data aa wfliH as teats have 
demonatratad the fenitafltty of iapOeaentlng an operatioaal alNOttOim for 
eatiaation of inaoilation on the baala of polarortiting aatePlte dhita. data 

or oth«* aatenite data oould be used as wi3L Cloud amount and Olotki type are the 
{siaary predtotora. Total waUr vapor la audi laas algilfloant statAittfiaTly, but the 
"water vapoc** term is aMlad to inctlute airfuse preaairo (or alternatively, hei^ 
above aea laveiD as weU as lOant path (through aeo Z). Whan dhmA variations are 
excluded in analysis, there is au^ leas variant in the water tern. To some 
extent there is also a bu£U>^ oompenaatlon on the optiMlthiokiMaBt the larger the 
average alant path thte to latibute, tiie saaHar is the oi^ioel ttdokneaa due to 
water vapor. 

Generally, tiie seleotim of parametera has not led to any obvl^ ooncluaion 
on the apatM bias, if any. Some exaninatian of r^eaaion ooefTLoLentB has shown 
no oc^MTont pattern of a geogri^Moal or alte-epeaifLo d^endenoy. On ottMr 
hand, any reoonaideratian of agiatial bias will requlra a doaer look at ttM 
geop'^^Moal variatlm of diumai aiwi seasonal insclstion distributiras, whloh oan 
d^art fhoB noraal distributions. A s e asona l sepa rat ion of r e g resa i oo ooeftloients 
has already been adopted. 

The p res ent teohnLqiM makes use of a set of regresBlni coefflcimts for each 
hour (dtslng any given season) to provide an estimate of the total daily insolation 
for the dosesb daytime period. For appUoation to any lo^tloi covered by a polar- 
orbiting satellite, aH that is req uir ed, in additioo to the pertinent houriy recp^esBlon 
ooeCOoimts, is the time, latitude, and solar declination (aasumit^ that the ii^ldent 
solar flux acroas a horizontal surface at the top of the atmosphere Is known). 
Individual houriy estimates are ooabined through wei^iting fkctors that depend on 
the existing and maxiBuB solar zeniih angles. 

The p res ent method require e that the analyst supply three pdeces of 
informations the gyoae doud type, the doud amount, and the totd predpttabLe 
water or aurface dewpoint. DniUally the doud type and amount will be determined 
by an analyst using an interaotive system displaying aataUlte imagery. The 
feeatbOity of thia approach waa temonstnted in the interactive tests of Section 
5.4. If not available thorn satellite sounder profUse of rnixix^ ratio, the 
predpitable water quantity oan be acquired thorn NMC analyses or thorn 
dimatdogioal values if no other source is available. For eadi time of satellite 
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oov««g«, tabulated ragnoBion ooafOoimtB art appliad to tha thraa paraattva to 
obtain a^fflatas of dBDy InBolation. Raaults for tha aavaral hom ara waigbtad to 
ani.va at a daily aatimata. Tha intvaotlva anOyab oan uaa his knowladga of 
atmoapharlo disturbanoaa and thatr tfiaraotaiatio dloud aignaturaa to altar 
wal^iting fotcra if ha baUavaa Uw avaUahla cdoud ^aervationa do not aoomtaly 
r a pr aaa nt Vtm day as a whcOa fhr any partioular location. Slailai^, tha talas 
introduoad tay ttw rag ^earton pr o om t whL^ results in tha owastimatiaa of 
insolation on mo^ doudy days and tha urktareatioatlon of insolation on u.oatly 
dear days oan ba ramovad if tha analyst dataots althar of thasa two dtaatiofa and 
oompanaatas acoordinipy. 

It is Ukaly that algorithm improvamants oan ba aooompUshad with improvad 
paramatailzation and weighting. As long as surfkoa data are used in ttM estimation 
prooaaa» surfaoe dawpotots oan ba uaad (as aatlmatarB of total pracdpitada water). 
However, an alternative is required ftom aateUita data. At tha p rese nt time it 
igpeera that the retrieved mixing ratioB- that wn indudad in tha TOYS product 
should be appliad, but other measures might be found also. Mors work is required 
on tiM doud parametera, with poaaible omiaLdtaratlm of an automated treatment of 
tha larger scale doud orgardzatim ChdpfULfor flUlng in gaps between pdar-orbiter 
data ooverage). Refleotance data, both of doud and of the siaDaoe, should be 
introduced even if limited to a slne^ time eadi day. When availade, the diurnal 
range of- temperature oan be used as an additimal estimator of dally Insdatim or 
at least as a ohaok on the estimated insolation. 

One step in the improvement prooaas is to engage in iterativa ragraaaim to a 
greater extent that was done here, m other words, first solutions of ooefELoients 
are used to infer doud parameters, aik! the infairad doud parameters are tlMn 
used to d^ain new coefOeients, with tha process continuing until a convergence la 
obtained. Ftalhermore, if the final results tend to be nohlineu*, departures of the 
estimates fhom the original observations can be used in another re gr e s s i on to 
diminate tiie bias. 

As tha pdar<rbiting satellite data are used in the algorithm, it wilL be 
beneficial in tiM learning stage to match the aatdlite data with tha SOLMET data. 
This wm provide the most dirwt link to th* oeslrad parameters. Whm InteractiLva 
procedures are used, it may be most profitable for the analyst to estimate 
parameters over the sites of interest at other times than those of satellite possaga. 
In fkct, the times can ba fhrad near vooa, or at critical d epa rtures fhom noon so 
that fixed weigits can be ^>plied. 
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OttMT ttUiDite data, at td^ raanilutlon, ooidd ba uaad to aotiaata a apaoiflo 
oiloud aaotait and taaparatira. Raaulta for mch grid liLaaant ooidd ba atorad with 
tha TOYS produot aua aary. &rfhoa raflmUinoa di^ ahoiild ba a<Mad to p«*aJt an 
luldttiKMdoorraotioo to tha inadlation aatiaation wImd naoaaaary (a^., with a fiaah 
anow oovar). 

Finally, tha SOLMET atationa witti both dlraot and difSVaa adlar radiation 
aaaaraaanta ahouLd ba aaainad aaparatdly to aradda a diraot-difXVaw inadianoa 
a^aration. Thia aa pa ra t l on Aould laad uMaataily to a battar inaoilation aatlaatlon 
algcritha, aapanlaTly fbr otloudy and partly doudy oonditiona. 

It ia raoonaaoctod that aotion ba takan to ooaidata tha <tovalopmant of tha 
uaa of tha HIRS data in tha apaoifloatian of cloud paramatare. Raaidta or 
auaaariaa dartvad tharafi?oa, avan if not uaad in an automatad analy^ would aid 
in tiM dadaion Baking pr ooeaa of tha intaraotlva andyat. 
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APPENDnt HIRS DATA ANALYSIS FOR CLOUD TYPE 


As a first st^ to the analysis of HIRS data into daases relatad to doud 
parametsrs of intarest for tiw estimation of insoilationy ciloud tarigata had to be 
setLsoted on the basis of known doud ^niMS. Prints GOES imagery were acquired 
to assist is the selactlm of the doud targets, with Initial emphasb ca doud ^pe 
rather than qpeoiflc amount. These seleotlons r e prese n t the first i^iase of a 
supervised duster analysis. In effect, ttie doud ig)eciflcati«)s, in terms of 
previoudy establiatod doud type categories, constitute d the dxxid truth for tte 
study. 

From all available HIRS infiered data, seven channels, centered at wave 
numbers of 705, 732, 7^, 901, 1218, 2191 and 2210 cm**‘l, were selected initially 
for flrOier analysis. An effort was made to avoid excessive duplication in channel 
respo ns e, but to maintain a minimum number diannela consirtent with the 
described number of doud types. Data fhom eadi selected channel were limb 
corrected; selected printouts were made in terms of scadixMS and scaiapots. The 
ap p licatlMi of latitude and Imgltude grid lines to the printouts makes it possible to 
dieck the accuracy data locations and matdnqs. A gray-scaled printout of 
infirared win(k>w data fhom HIRS can be diecked readily against the gridded image 
of the GOES infrared window. 

Diitial doud types were defined so as to avoid ambiguities as much as 
possible. Six adopted types were CLEAR, HIGH, MID, LOW, MIXED and THICK. A 
stepwise .muMple discriminant analysis was sq>plied to the HIRS data (uring the 
sevm diannds as predictors to specify six doud types). Table A-1 sum marizes the 
results obtained fbom the SPSS discriminant analysis package as ;q)plied to a June 
26,1980 data swath(N0AA-6 satdlite). It is apparent that dear sides over codL 
water were interpreted as low douds in a number at cases. The HIGH category 
shows same apraai associated with variations in opacity. As mi^ be expected, 
the MIXED category sometimes is assigned to a particular cloud type. 

Figure A-1 is a plot in the ^ce of the first two discriminant flioctlons. Sets 
of coefficients derived for each fbnctlcxi are applied to each observaticxi in 
(teterminatlcxi of the category. Separation of data points, each assigned a doud 
category, Is eq>parmit la Figure A-1, but distances betwem groups are not large. 
The same discriminant fbnctlcxa were applied in an independent test to HIRS data 
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TABLE A-1 


SUMMARY OF INITIAL CLOUD CLASSIFICATIONS FROM 
HIRS/2 MULTISPECTRALSCANSPOT DATA 


CLASSEFICAHON RESULTS 


ACTUAL GROUP 

NO. OF 
CASES 

PREDICTED GROUP MEMBERSHIP 
12 3 4 

5 

6 

CLEAR 

1 

62 

51 

0 

0 

0 

0 

0 




82.3 

0.0 

0.0 

17.7 

0.0 

0.0 

HIGH 

*2 

24 

1 

19 

2 

2 

0 

0 




4.2 

79.2 

8.3 

8.3 

0.0 

0.0 

MID 

3 

17 

0 

0 

16 

0 

1 

0 




0.0 

0.0 

94.1 

0.0 

5.9 

0.0 

^OW 

4 

22 

0 

0 

0 

22 

0 

0 




0.0 

0.0 

0.0 

100.0 

0.0 

0.0 

MIXED 

5 

43 

0 

1 

2 

1 

39 

0 




0.0 

2.3 

4.7 

2.3 

90.7 . 

0.0 

THICK 

6 

32 

0 

0 

0 

0 

1 

31 




0.0 

0.0 

0.0 

0.0 

3.1 

96.9 


PERCENT CF GROUPED CASES CORRECTLY CLASSIFIED 89.00 
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FIGURE A-1 CLOUD CLASSIFICATION DISTRIBUTION BASED ON 
HIRS DATA FOR FIRST TWO DISCRIMINANT FUNCTIONS , 



firoB 00^>-0215 GMT, 30 June I960, ow aiatarn North ABodoi and the coastal 
AUantio* A ooBp30te diatribution of predioted doud ^rpea is atown in the printout 
diqilayed in Fig. A-2. There is UtOe hi^s-qualtty ground truth avaGahle tor 
dieoldng the predioted ^tialdiatrtbutloa of types. Some spot dieoks are poastbiLe 
fhom NortlMm HeoisplMrio surface charts (not reproduced here). Several such 
reports, althoutfi not strictly coinoidenb with HIRS data, substantiated significant 
dlasBifloations. Figire A-3 is a copy of a viatble print from GOES imagery ft^m 
2315Z, 29 June I960, approximately one and one-half hom prior to the satellite 
swath over VbB Eastern United States (0030 to 0045Z, 30 June 80). Later satellite 
imagery applies after sunset over the Eastern portion of the area. Several features 
apparttit in the satellite imagery are well represttited in Figure A-2. The large 
area of thunderstorms over northern Miasiasippi ard Alabama is evident as a cluster 
of thick clouds (category 6) in Figure A-2. The predpitatiai producing clouds over 
New York extending into New England are deariy described. The large swiil of 
Boatlylow doudineas over the Great Lakes is well represented as is the large dear 
area in the Ohio Valley. The dear area or "dry tongue" extending fhom Western 
New York aouthweabward into Kentudcy can be seen readily. Althoufd^ other 
fsatures are difficult to ddacem in the iUustreted copy of the GOES image, a 
qualitative evaluation of the dcud type estimatiai modd. using the original print 
indioates c(X)sLderabla success. No doubt improvement would result ftcm the use 
of additional multispectral parametans. Therefore, this method should be fdly 
developed for application. 
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FIGURE A- 2 ESTIMATED CLOUD CATEGORIES FROM INDEPENDENT 
HIRS SCANSPOTS. 29 JUNE 1980 
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FIGURE A-3 COPY OF GOES VISIBLE IMAGE NEAR SUNSET. 29 JUNE 1980 
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